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ABSTRACT: The physical conditions in molecular clouds control the nature and rate of 
star formation, with consequences for planet formation and galaxy evolution. The focus of 
this review is on the conditions that characterize regions of star formation in our Galaxy. A 
review of the tools and tracers for probing physical conditions includes summaries of generally 
applicable results. Further discussion distinguishes between the formation of low-mass stars 
in relative isolation and formation in a clustered environment. Evolutionary scenarios and 
theoretical predictions are more developed for isolated star formation, and observational tests are 
beginning to interact strongly with the theory. Observers have identified dense cores collapsing 
to form individual stars or binaries, and analysis of some of these support theoretical models of 
collapse. Stars of both low and high mass form in clustered environments, but massive stars form 
almost exclusively in clusters. The theoretical understanding of such regions is considerably less 
developed, but observations are providing the ground rules within which theory must operate. 
The most rich and massive star clusters form in massive, dense, turbulent cores, which provide 
models for star formation in other galaxies. 
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1 INTRODUCTION 



Long after their parent spiral galaxies have formed, stars continue to form by 
repeated condensation from the interstellar medium. In the process, parts of 
the interstellar medium pass through a cool, relatively dense phase with a great 
deal of complexity- molecular clouds. While both the diffuse interstellar medium 
and stars can be supported by thermal pressure, most molecular clouds cannot 
be thermally supported (Goldreich & Kwan 1974). Simple consideration would 
suggest that molecular clouds would be a very transient phase in the conversion 
of diffuse gas to stars, but in fact they persist much longer than expected. During 
this extended life, they produce an intricate physical and chemical system that 
provides the substrate for the formation of planets and life, as well as stars. 
Comparison of cloud masses to the total mass of stars that they produce indicates 
that most of the matter in a molecular cloud is sterile; stars form only in a small 
fraction of the mass of the cloud (Leisawitz et al 1989). 

The physical conditions in the bulk of a molecular cloud provide the key to 
understanding why molecular clouds form an essentially metastable state along 
the path from diffuse gas to stars. Most of the mass of most molecular clouds 
in our Galaxy is contained in regions of modest extinction, allowing photons 
from the interstellar radiation field to maintain sufficient ionization for magnetic 
fields to resist collapse (McKee 1989); most of the molecular gas is in fact in 
a photon-dominated region, or PDR (Hollenbach & Tielens 1997). In addition, 
most molecular gas has supersonic turbulence (Zuckerman & Evans 1974). The 
persistence of such turbulence over the inferred lifetimes of clouds in the face of 
rapid damping mechanisms (Goldreich & Kwan 1974) suggests constant replen- 
ishment, most likely in a process of self-regulated star formation (Norman & Silk 
1980, Bertoldi & McKee 1996), since star formation is accompanied by energetic 
outflows, jets, and winds (Bachiller 1996). 

For this review, the focus will be on the physical conditions in regions that are 
forming stars and likely precursors of such regions. While gravitational collapse 
explains the formation of stars, the details of how it happens depend critically 
on the physical conditions in the star-forming region. The details determine the 
mass of the resulting star and the amount of mass that winds up in a disk around 
the star, in turn controlling the possibilities for planet formation. The physical 
conditions also control the chemical conditions. With the recognition that much 
interstellar chemistry is preserved in comets (Crovisier 1999, van Dishoeck & 
Blake 1998), and that interstellar chemistry may also affect planet formation and 
the possibilities for life (e.g. Pendleton 1997, Pendleton & Tielens 1997, Chyba 
&: Sagan 1992), the knowledge of physical conditions in star-forming regions has 
taken on additional significance. 
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Thinking more globally, different physical conditions in different regions deter- 
mine whether a few, lightly clustered stars form (the isolated mode) or a tight 
grouping of stars form (the clustered mode) (Lada 1992; Lada et al 1993). The 
star formation rates per unit mass of molecular gas vary by a factor > 10 2 in 
clouds within our own Galaxy (Evans 1991, Mead et al 1990), and starburst 
galaxies achieve even higher rates than are seen anywhere in our Galaxy (e.g. 
Sanders et al 1991). Ultimately, a description of galaxy formation must incor- 
porate an understanding of how star formation depends on physical conditions, 
gleaned from careful study of our Galaxy and nearby galaxies. 

Within the space limitations of this review, it is not possible to address all the 
issues raised by the preceding overview. I will generally avoid topics that have 
been recently reviewed, such as circumstellar disks (Sargent 1996, Papaloizou & 
Lin 1995, Lin & Papaloizou 1996, Bodenheimer 1995), as well as bipolar out- 
flows (Bachiller 1996), and dense PDRs (Hollenbach &: Tielens 1997). I will also 
discuss the sterile parts of molecular clouds only as relevant to the process that 
leads some parts of the cloud to be suitable for star formation. While the chem- 
istry and physics of star-forming regions are coupled, chemistry has been recently 
reviewed (van Dishoeck & Blake 1998). Astronomical masers are being concur- 
rently reviewed (Menten 1999); as with HII regions, they will be discussed only 
as signposts for regions of star formation. 

I will focus on star formation in our Galaxy. Nearby regions of isolated, low- 
mass star formation will receive considerable attention (§|j) because we have 
made the most progress in studying them. Their conditions will be compared 
to those in regions forming clusters of stars, including massive stars (§||). These 
regions of clustered star formation are poorly understood, but they probably form 
the majority of stars in our Galaxy (Elmegreen 1985), and they are the regions 
relevant for comparisons to other galaxies. 

Even with such a restricted topic, the literature is vast. I make no attempt 
at completeness in referencing. On relatively non-controversial topics, I will tend 
to give an early reference and a recent review; for more unsettled topics, more 
references, with different points of view, will be given. Recent or upcoming 
publications with significant overlap include Hartmann (1998), Lada & Kylafis 
(1999), and Mannings et al (2000). 

2 PHYSICAL CONDITIONS 

The motivation for studying physical conditions can be found in a few simple 
theoretical considerations. Our goal is to know when and how molecular gas 
collapses to form stars. In the simplest situation, a cloud with only thermal 
support, collapse should occur if the mass exceeds the Jeans (1928) mass, 

Mj = L V - 5 = 18M T^- 5 n- - 5 , (1) 

where Tk is the kinetic temperature (kelvins), p is the mass density (gm cm -3 ), 
and n is the total particle density (cm~ 3 ). In a molecular cloud, H nuclei are 
almost exclusively in H2 molecules, and n ~ n(H2) + n(He). Then p = p, n mHn, 
where run is the mass of a hydrogen atom and p n is the mean mass per particle 
(2.29 in a fully molecular cloud with 25% by mass helium). Discrepancies between 
coefficients in the equations presented here and those in other references usually 
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are traceable to a different definition of n. In the absence of pressure support, 
collapse will occur in a free-fall time (Spitzer 1978), 

/ 3tt \ - 5 

tff = i^G-p) = 3 - 4 x lo7 ^ a5 y ears - ( 2 ) 

If Tk = 10 K and n > 50 cm -3 , typical conditions in the sterile regions (e.g. 
Blitz 1993), Mj < 80M Q , and tff < 5 x 10 6 years. Our Galaxy contains about 
1-3 x 10 9 M of molecular gas (Bronfman et al 1988, Clemens et al 1988, Combes 
1991). The majority of this gas is probably contained in clouds with M > 10 4 
M Q (Elmegreen 1985). It would be highly unstable on these grounds, and free- 
fall collapse would lead to a star formation rate, M* > 200 M yr _1 , far in 
excess of the recent Galactic average of 3 M yr _1 (Scalo 1986). This argument, 
first made by Zuckerman & Palmer (1974), shows that most clouds cannot be 
collapsing at free fall (see also Zuckerman & Evans 1974). Together with evidence 
of cloud lifetimes of about 4 x 10 7 yr (Bash et al 1977; Leisawitz et al 1989), this 
discrepancy motivates an examination of other support mechanisms. 

Two possibilities have been considered, magnetic fields and turbulence. Calcu- 
lations of the stability of magnetized clouds (Mestel & Spitzer 1956, Mestel 1965) 
led to the concept of a magnetic critical mass (Mb). For highly flattened clouds, 

M B = (2vr)- 1 G-°- 5 $, (3) 
(Li k, Shu 1996), where $ is the magnetic flux threading the cloud, 



<I> 



= j Bda. (4) 



Numerical calculations (Mouschovias & Spitzer 1976) indicate a similar coefficient 
(0.13). If turbulence can be thought of as causing pressure, it may be able to 
stabilize clouds on large scales (e.g. Bonazzola et al 1987). It is not at all clear 
that turbulence can be treated so simply. In both cases, the cloud can only be 
metastable. Gas can move along field lines and ambipolar diffusion will allow 
neutral gas to move across field lines with a timescale of (McKee et al 1993) 

3 

tAD = -r-^ ^ 7.3 x 10 13 x e years, (5) 

where r n j is the ion-neutral collision time. The ionization fraction (x e ) depends 
on ionization by photons and cosmic rays, balanced by recombination. It thus 
depends on the abundances of other species (X(x) = n(x)/n). 

These two suggested mechanisms of cloud support (magnetic and turbulent) 
are not entirely compatible because turbulence should tangle the magnetic field 
(compare the reviews by Mouschovias 1999 and McKee 1999). A happy mar- 
riage between magnetic fields and turbulence was long hoped for; Arons and 
Max (1975) suggested that magnetic fields would slow the decay of turbulence 
if the turbulence was sub-Alfvenic. Simulations of MHD turbulence in systems 
with high degrees of symmetry supported this suggestion and indicated that 
the pressure from magnetic waves could stabilize clouds (Gammie & Ostriker 
1996). However, more recent 3-D simulations indicate that MHD turbulence de- 
cays rapidly and that replenishment is still needed (Mac Low et al 1998, Stone et 
al 1999). The usual suggestion is that outflows generate turbulence, but Zweibel 
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(1998) has suggested that an instability induced by ambipolar diffusion may con- 
vert magnetic energy into turbulence. Finally, the issue of supporting clouds 
assumes a certain stability and cloud integrity that may be misleading in a dy- 
namic interstellar medium (e.g. Ballesteros-Paredes et al 1999). For a current 
review of this field, see Vazquez- Semadini et al (2000). 

With this brief and simplistic review of the issues of cloud stability and evolu- 
tion, we have motivated the study of the basic physical conditions: 

T K ,n,v,B,X. (6) 

All these are local variables; in principle, they can have different values at each 
point (r) in space and can vary with time (i). In practice, we usually can mea- 
sure only one component of vector quantities, integrated through the cloud. For 
example, we measure only line-of-sight velocities, usually characterized by the 
linewidth (Av) or higher moments, and the line-of-sight magnetic field (B z ) 
through the Zeeman effect, or the projected direction (but not the strength) 
of the field in the plane of the sky (B±) by polarization studies. In addition, 
our observations always average over finite regions, so we attempt to simplify 
the dependence on r by assumptions or models of cloud structure. In §^, I will 
describe the methods used to probe these quantities and some overall results. 
Abundances have been reviewed recently (van Dishoeck & Blake 1998), so only 
relevant results will be mentioned, most notably the ionization fraction, x e . 

In addition to the local variables, quantities that explicitly integrate over one 
or more dimensions are often measured. Foremost is the column density, 

N = f ndl. (7) 



The extinction or optical depth of dust at some wavelength is a common surrogate 
measure of N . If the column density is integrated over an area, one measure of 
the cloud mass within that area is obtained: 

M N = J Nda. (8) 

Another commonly used measure of mass is obtained by simplification of the 
virial theorem. If external pressure and magnetic fields are ignored, 

My = C v G~ x RAv 2 = 2WM Q C v R{pc){Av(km s" 1 )) 2 , (9) 

where R is the radius of the region, Av is the FWHM linewidth, and the constant 
(C v ) depends on geometry and cloud structure, but is of order unity (e.g. McKee 
& Zweibel 1992, Bertoldi & McKee 1992). A third mass estimate can be obtained 
by integrating the density over the volume, 

M n = [ ndv. (10) 



M n is commonly used to estimate f v , the volume filling factor of gas at density n, 
by dividing M n by another mass estimate, typically My. Of the three methods 
of mass determination, the virial mass is the least sensitive to uncertainties in 
distance and size, but care must be taken to exclude unbound motions, such as 
outflows. 

Parallel to the physical conditions in the gas, the dust can be characterized by 
a set of conditions: 
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T D ,n D ,K(u), 



(11) 



where Tp is the dust temperature, no is the density of dust grains, and k(u) 
is the opacity at a given frequency, v. If we look in more detail, grains have a 
range of sizes (Mathis et al 1977) and compositions. For smaller grains, Td is 
a function of grain size. Thus, we would have to characterize the temperature 
distribution as a function of size, the composition of grains, both core and mantle, 
and many more optical constants to capture the full range of grain properties. 
For our purposes, To and k(u) are the most important properties, because Tp 
affects gas energetics and To and k(u) control the observed continuum emission 
of molecular clouds. The detailed nature of the dust grains may come into play 
in several situations, but the primary observational manifestation of the dust is 
its ability to absorb and emit radiation. For this review, a host of details can be 
ignored. The optical depth is set by 



where it is convenient to define k(v) so that N is the gas column density rather 
than the dust column density. Away from resonances, the opacity is usually 
approximated by k{i>) oc u 13 . 



The most fundamental fact about molecular clouds is that most of their contents 
are invisible. Neither the H2 nor the He in the bulk of the clouds are excited 
sufficiently to emit. While fluorescent emission from H2 can be mapped over 
the face of clouds (Luhman & Jaffe 1996), the ultraviolet radiation needed to 
excite this emission does not penetrate the bulk of the cloud. In shocked regions, 
H2 emits rovibrational lines that are useful probes of T% and v(r) (e.g. Draine 
& McKee 1993). Absorption by H2 of background stars is also difficult: the 
dust in molecular clouds obnubilates the ultraviolet that would reveal electronic 
transitions; the rotational transitions are so weak that only huge N would pro- 
duce absorption, and the dust again obnubilates background sources; only the 
vibrational transitions in the near-infrared have been seen in absorption in only 
a few molecular clouds (Lacy et al 1994). Gamma rays resulting from cosmic 
ray interactions with atomic nuclei do probe all the material in molecular clouds 
(Bloemen 1989, Strong et al 1994). So far, gamma-ray studies have suffered from 
low spatial resolution and uncertainties in the cosmic ray flux; they have been 
used mostly to check consistency with other tracers on large scales. 

In the following subsections, I will discuss probes of different physical quanti- 
ties, including some general results, concluding with a discussion of the observa- 
tional and analytical tools. Genzel (1992) has presented a detailed discussion of 
probes of physical conditions. 

3.1 Tracers of Column Density, Size, and Mass 

Given the reticence of the bulk of the H2, essentially all probes of physical condi- 
tions rely on trace constituents, such as dust particles and molecules other than 
H2. Dust particles (Mathis 1990, Pendleton & Tielens 1997) attenuate light at 
short wavelengths (ultraviolet to near-infrared) and emit at longer wavelengths 



t d (u) = k(u)N, 



(12) 
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(far-infrared to millimeter) . Assuming that the ratio of dust extinction at a fixed 
wavelength to gas column density is constant, one can use extinction to map 
N in molecular clouds, and early work at visible wavelengths revealed locations 
and sizes of many molecular clouds before they were known to contain molecules 
(Barnard 1927, Bok Sz Reilly 1947, Lynds 1962). There have been more recent 
surveys for small clouds (Clemens & Barvainis 1988) and for clouds at high lat- 
itude (Blitz et al 1984). More recently, near-infrared surveys have been used to 
probe much more deeply; in particular, the H — K color excess can trace N to an 
equivalent visual extinction, Ay ~ 30 mag (Lada et al 1994, Alves et al 1998a). 
This method provides many pencil beam measurements through a cloud toward 
background stars. The very high resolution, but very undersampled, data require 
careful analysis but can reveal information on mass, large-scale structure in iV 
and unresolved structure (Alves et al 1998a). Padoan et al (1997) interpret the 
data of Lada et al (1994) in terms of a lognormal distribution of density, but 
Lada et al (1999) show that a cylinder with a density gradient, n(r) cc r~ 2 , also 
matches the observations. 

Continuum emission from dust at long wavelengths is complementary to ab- 
sorption studies (e.g. Chandler & Sargent 1997). Because the dust opacity de- 
creases with increasing wavelength (k(u) cx u 13 , with (3 ~ 1 — 2), emission at 
long wavelengths can trace large column densities and provide independent mass 
estimates (Hildebrand 1983). The data can be fully sampled and have reason- 
ably high resolution. The dust emission depends on the dust temperature (To), 
linearly if the observations are firmly in the Rayleigh- Jeans limit, but exponen- 
tially on the Wien side of the blackbody curve. Observations on both sides of the 
emission peak can constrain Try. Opacities have been calculated for a variety of 
scenarios including grain mantle formation and collisional concrescence (e.g. Os- 
senkopf & Henning 1994). For grain sizes much less than the wavelength, [3 ~ 2 
is expected from simple grain models, but observations of dense regions often 
indicate lower values. By observing at a sufficiently long wavelength, one can 
trace N to very high values. Recent results indicate that ro(1.2mm) = 1 only 
for Ay ~ 4 x 10 4 mag (Kramer et al 1998a) in the less dense regions of molec- 
ular clouds. In dense regions, there is considerable evidence for increased grain 
opacity at long wavelengths (Zhou et al 1990, van der Tak et al 1999), suggesting 
grain growth through collisional concrescence in addition to the formation of icy 
mantles. Further growth of grains in disks is also likely (Chandler h Sargent 
1997). 

The other choice is to use a trace constituent of the gas, typically molecules that 
emit in their rotational transitions at millimeter or submillimeter wavelengths. 
By using the appropriate transitions of the appropriate molecule, one can tune 
the probe to study the physical quantity of interest and the target region along 
the line of sight. This technique was first used with OH (Barrett et al 1964), and 
it has been pursued intensively in the 30 years since the discovery of polyatomic 
interstellar molecules (Cheung et al 1968; van Dishoeck & Blake 1998). 

The most abundant molecule after H2 is carbon monoxide; the main isotopomer 
( 12 C 16 0) is usually written simply as CO. It is the most common tracer of molec- 
ular gas. On the largest scales, CO correlates well with the gamma-ray data 
(Strong et al 1994), suggesting that the overall mass of a cloud can be measured 
even when the line is quite opaque. This stroke of good fortune can be under- 
stood if the clouds are clumpy and macroturbulent, with little radiative coupling 
between clumps (Wolfire et al 1993); in this case, the CO luminosity is propor- 
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tional to the number of clumps, hence total mass. Most of the mass estimates 
for the larger clouds and for the total molecular mass in the Galaxy and in other 
galaxies are in fact based on CO. On smaller scales, and in regions of high column 
density, CO fails to trace column density, and progressively rarer isotopomers are 
used to trace progressively higher values of N. Dickman (1978) established a 
strong correlation of visual extinction Ay with 13 CO emission for 1.5 < Ay < 5. 
Subsequent studies have used C 18 and C 17 to trace still higher N (Frerking 
et al 1982). These rarer isotopomers will not trace the outer parts of the cloud, 
where photodissociation affects them more strongly than the common ones, but 
we are concerned with the more opaque regions in this review. 

When comparing N measured by dust emission with N traced by CO iso- 
topomers, it is important to correct for the fact that emission from low-J transi- 
tions of optically thin isotopomers of CO decreases with Tk, while dust emission 
increases with Td (Jaffe et al 1984). Observations of many transitions can avoid 
this problem but are rarely done. To convert iV(CO) to N requires knowledge 
of the abundance, X(CO). The only direct measure gave X(CO) = 2.7 x 10~ 4 
(Lacy et al 1994), three times greater than inferred from indirect means (e.g. 
Frerking et al 1982). Clearly, this area needs increased attention, but at least 
a factor of three uncertainty must be admitted. Studies of some particularly 
opaque regions in molecular clouds indicate severe depletion (Kuiper et al 1996, 
Bergin & Langer 1997), raising the concern that even the rare CO isotopomers 
may fail to trace N. Indeed, Alves et al (1999) find that C 18 fails to trace 
column density above Ay = 10 in some regions, and Kramer et al (1999) argue 
that this failure is best explained by depletion of C 18 0. 

Sizes of clouds, characterized by either a radius (R) or diameter (I), are mea- 
sured by mapping the cloud in a particular tracer; for non-spherical clouds, these 
are often the geometric mean of two dimensions, and the aspect ratio (a/b) char- 
acterizes the ratio of long and short axes. The size along the line of sight (depth) 
can only be constrained by making geometrical assumptions (usually of a spheri- 
cal or cylindrical cloud). One possible probe of the depth is H3 + , which is unusual 
in having a calculable, constant density in molecular clouds. Thus, a measure- 
ment of A^Hg") can yield a measure of cloud depth (Geballe & Oka 1996). 

With a measure of size and a measure of column density, the mass (Mat) may 
be estimated (equation 8); with a size and a linewidth (Av), the virial mass 
(My) can be estimated (equation 9). On the largest scales, the mass is often 
estimated from integrating the CO emission over the cloud and using an empirical 
relation between mass and the CO luminosity, L(CO). The mass distribution 
has been estimated for both clouds and clumps within clouds, primarily from 
CO, 13 CO, or C 18 0, using a variety of techniques to define clumps and estimate 
masses (Blitz 1993, Kramer et al 1998b, Heyer Sz Terebey 1998, Heithausen et 
al 1998). These studies have covered a wide range of masses, with Kramer et al 
extending the range down to M = 1O -4 M0. The result is fairly well agreed on: 
dN(M) ex M- a dM, with 1.5 < a < 1.9. Elmegreen & Falgarone (1996) have 
argued that the mass spectrum is a result of the fractal nature of the interstellar 
gas, with a fractal dimension, D = 2.3 ±0.3. There is disagreement over whether 
clouds are truly fractal or have a preferred scale (Blitz & Williams 1997). The 
latter authors suggest a scale of 0.25-0.5 pc in Taurus based on 13 CO. On the 
other hand, Falgarone et al (1998), analyzing an extensive data set, find evidence 
for continued structure down to 200 AU in gas that is not forming stars. The 
initial mass function (IMF) of stars is steeper than the cloud mass distribution 
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for M+ > 1M© but is natter than the cloud mass function for M* < 1M Q (e.g. 
Scalo 1998). Understanding the origin of the differences is a major issue (see 
Williams et al 2000, Meyer et al 2000). 

3.2 Probes of Temperature and Density 

The abundances of other molecules are so poorly constrained that CO isotopomers 
and dust are used almost exclusively to constrain N and M^. Of what use are 
the over 100 other molecules? While many are of interest only for astrochemistry, 
some are very useful probes of physical conditions like Tk,u,v,B z , and x e . 

Density (n) and gas temperature (Tk) are both measured by determining the 
populations in molecular energy levels and comparing the results to calculations 
of molecular excitation. A useful concept is the excitation temperature (T ex ) 
of a pair of levels, defined to be the temperature that gives the actual ratio of 
populations in those levels, when substituted into the Boltzmann equation. In 
general, collisions and radiative processes compete to establish level populations; 
when lines are optically thick, trapping of line photons enhances the effects of 
collisions. For some levels in some molecules, radiative rates are unusually low, 
collisions dominate, T cx = Tk, and observational determination of these "ther- 
malized" level populations yields Tk- Un-thermalized level populations depend 
on both n and Tk', with a knowledge of Tk, observational determination of these 
populations yields n, though trapping usually needs to be accounted for. While 
molecular excitation probes the local n and Tk in principle, the observations 
themselves always involve some average over the finite beam and along the line 
of sight. Consequently, a model of the cloud is needed to interpret the obser- 
vations. The simplest model is of course a homogeneous cloud, and most early 
work adopted this model, either explicitly or implicitly. 

Tracers of temperature include CO, with its unusually low dipole moment, and 
molecules in which transitions between certain levels are forbidden by selection 
rules. The latter include different K ladders of symmetric tops like NH3, CH3CN, 
etc. (Ho & Townes 1983, Loren & Mundy 1984). Different ladders in H 2 CO 
also probe Tk in dense, warm regions (Mangum and Wootten 1993). A useful 
feature of CO is that its low- J transitions are both opaque and thermalized 
in most parts of molecular clouds. In this case, observations of a single line 
provide the temperature, after correction for the cosmic background radiation 
and departures from the Rayleigh-Jeans approximation (Penzias et al 1972, Evans 
1980). Early work on CO (Dickman 1975) and NH 3 (Martin & Barrett 1978) 
established that Tk — 10 K far from regions of star formation and that sites of 
massive star formation are marked by elevated Tk, revealed by peaks in maps of 
CO (e.g. Blair et al 1975). 

The value of Tk far from local heating sources can be understood by balanc- 
ing cosmic ray heating and molecular cooling (Goldsmith & Langer 1978), while 
elevated values of Tk in star forming regions have a more intricate explanation. 
Stellar photons, even when degraded to the infrared, do not couple well to molec- 
ular gas, so the heating goes via the dust. The dust is heated by photons and 
the gas is heated by collisions with the dust (Goldreich & Kwan 1974); above 
a density of about 10 4 cm' 3 , Tk becomes well coupled to To (Takahashi et 
al 1983). Observational comparison of Tk to Tp, determined from far- infrared 
observations, supports this picture (e.g. Evans et al 1977, Wu & Evans 1989). 

In regions where photons in the range of 6 to 13.6 eV impinge directly on 
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molecular material, photoelectrons ejected from dust grains can heat the gas 
very effectively and Tk may exceed Tjj. These PDRs (Hollenbach & Tielens 
1997) form the surfaces of all clouds, but the regions affected by these photons 
are limited by dust extinction to about Ay ~ 8 mag (McKee 1989). However, 
the CO lines often do form in the PDR regions, raising the question of why they 
indicate that Tk ~ 10K. Wolfire et al (1993) explain that the optical depth in 
the lower J levels usually observed reaches unity at a place where the Tk and 
n combine to produce an excitation temperature (T ex ) of about 10 K. Thus, the 
agreement of Tk derived from CO with the predictions of energetics calculations 
for cosmic ray heating may be fortuitous. The Tk derived from NH3 refer to 
more opaque regions and are more relevant to cosmic ray heating. Finally, in 
localized regions, shocks can heat the gas to very high Tk', values of 2000 K are 
observed in H2 ro- vibrational emission lines (Beckwith et al 1978). It is clear that 
characterizing clouds by a single Tk, which is often done for simplicity, obscures 
a great deal of complexity. 

Density determination requires observations of several transitions that are not 
in local thermodynamic equilibrium (LTE). Then the ratio of populations, or 
equivalently T ex , can be used to constrain density. A useful concept is the critical 
density for a transition from level j to level k, 

n c (jk) = A jk /j jk , (13) 

where Aj k is the Einstein A coefficient and njjk is the collisional deexcitation 
rate per molecule in level j. In general, both H2 and He are effective collision 
partners, with comparable collision rates, so that excitation techniques measure 
the total density of collision partners, n ~ n(H.2) +n(He). In some regions of high 
x e , collisions with electrons may also be significant. Detection of a particular 
transition is often taken to imply that n > n c {jk), but this statement is too 
simplistic. Lines can be seen over a wide range of n, depending on observational 
sensitivity, the frequency of the line, and the optical depth (e.g. Evans 1989). 
Observing high frequency transitions, multilevel excitation effects, and trapping 
all tend to lower the effective density needed to detect a line. Table 1 contains 
information for some commonly observed lines, including the frequency, energy 
in K above the effective ground state (E up (K)), and the critical densities at 
Tk = 10 K and 100 K. For comparison, the Table also has n e ff, the density 
needed to produce a line of 1 K, easily observable in most cases. The values 
of n e ff were calculated with a Large Velocity Gradient (LVG) code (§ |3.5| ) to 
account for trapping, assuming log(A r /Af ) = 13.5 for all species but NH3, for 
which log(N/Av) = 15 was used. N/Av has units of cm -2 (km s -1 ) -1 . These 
column densities are typical and produce modest optical depths. Note that n e ff 
can be as much as a factor of 1000 less than the critical density, especially for 
high excitation lines and high Tk- Clearly, the critical density should be used as 
a guideline only; more sophisticated analysis is necessary to infer densities. 

Assuming knowledge of Tk, at least two transitions with different n c {jk) are 
needed to determine both n and the line optical depth, r,-^ oc N k /Av, which 
determines the amount of trapping, and more transitions are desirable. Because 
Ajj-i oc J 3 , where J is the quantum number for total angular momentum, 
observing many transitions up a rotational energy ladder provides a wide range 
of n c {jk). Linear molecules, like HCN and HCO + , have been used in this way, 
but higher levels often occur at wavelengths with poor atmospheric transmission. 
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Relatively heavy species, like CS, have many accessible transitions, and up to 
five transitions ranging up to J = 10 have been used to constrain density (e.g. 
Carr et al 1995, van der Tak et al 1999). More complex species provide more 
accessible energy levels; transitions within a single K~i ladder of H2CO provide 
a valuable density probe (Mangum & Wootten 1993). Transitions of H2CO with 
AJ = are accessible to large arrays operating at centimeter wavelengths (e.g. 
Evans et al 1987). The lowest few of these H2CO transitions have the interesting 
property of absorbing the cosmic background radiation (Palmer et al 1969), Tex 
being cooled by collisional pumping (Townes & Cheung 1969). 

Application of these techniques to the homogeneous cloud model generally 
produces estimates of density exceeding 10 4 cm -3 in regions forming stars, while 
the sterile regions of the cloud are thought to have typical n ~ 10 2 — 10 3 cm -3 , 
though these are less well constrained. Theoretical simulations of turbulence 
have predicted lognormal (Vazquez-Semadini 1994) or power-law (Scalo et al 
1998) probability density functions. Studies of multiple transitions of different 
molecules with a wide range of critical densities often reveal evidence for density 
inhomogeneities; in particular, pairs of transitions with higher critical densities 
tend to indicate higher densities (e.g. Evans 1980, Plume et al 1997). Both 
density gradients and clumpy structure have been invoked to explain these results 
(see §H] and §|| for detailed discussion). Since lines with high n c (jk) are excited 
primarily at higher n, one can avoid to some extent the averaging over the line 
of sight by tuning the probe. 

3. 3 Kinematics 

In principle, information on v{r) is contained in maps of the line profile over the 
cloud. In practice, this message has been difficult to decode. Only motions along 
the line of sight produce Doppler shifts, and the line profiles average over the 
beam and along the line of sight. Maps of the line center velocity generally in- 
dicate that the typical cloud is experiencing neither overall collapse (Zuckerman 
& Evans 1974) nor rapid rotation (Arquilla & Goldsmith 1986, Goodman et al 
1993). Instead, most clouds appear to have velocity fields dominated by turbu- 
lence, because the linewidths are usually much greater than expected from ther- 
mal broadening. While such turbulence can explain the breadth of the lines, the 
line profile is not easily matched. Even a homogeneous cloud will tend to develop 
an excitation gradient in unthermalized lines because of trapping, and gradients 
in Tk or n toward embedded sources should exacerbate this tendency. Simple mi- 
croturbulent models with decreasing T ex (r) predict that self-reversed line profiles 
should be seen more commonly than they are. Models with many small clumps 
and macroturbulence have had some success in avoiding self-reversed line profiles 
(Martin et al 1984, Wolfire et al 1993, Falgarone et al 1994, Park & Hong 1995). 

The average linewidths of clouds are larger for larger clouds, the linewidth-size 
relation: Av oc ffl (Larson 1981). For clouds with the same N, the virial theorem 
would predict 7 = 0.5, consistent with the results of many studies of clouds as a 
whole (e.g. Solomon et al 1987). Myers (1985) summarized the different relations 
and distinguished between those comparing clouds as a whole and those studying 
trends within a single cloud. The status of linewidth-size relations within clouds, 
particularly in star-forming regions, will be discussed in later sections. 
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3.4 Magnetic Field and Ionization 

The magnetic field strength and direction are important but difficult to measure. 
Heiles et al (1993) review the observations and McKee et al (1993) review the- 
oretical issues. The only useful measure of the strength is the Zeeman effect, 
which probes the line-of-sight field, B z . Observations of HI can provide some 
useful probes of B z in PDRs (e.g. Brogan et al 1999), but cannot probe the 
bulk of molecular gas. Molecules suitable for Zeeman effect measurements have 
unpaired electrons and their resulting reactivity tends to decrease their abun- 
dance in the denser regions (e.g. Sternberg et al 1997). Almost all work has 
been done with OH, along with some work with CN (Crutcher et al 1996), but 
future prospects include CCS and excited states of OH and CH (Crutcher 1998). 
Measurements of B z have been made with thermal emission or absorption by OH 
(e.g. Crutcher et al 1993), mostly probing regions with n ~ 10 3 cm~ 3 , where 
B z ~ 20 uG or with OH maser emission, probing much denser gas, but with less 
certain conditions. As reviewed by Crutcher (1999a), the results for 14 clouds of 
widely varying mass with good Zeeman detections indicate that Mb is usually 
within a factor of 2 of the cloud mass. Given uncertainties, this result suggests 
that clouds with measured B z lie close to the critical-subcritical boundary (Shu 
et al 1999). The observations can be fit with B oc n , remarkably consistent 
with predictions of ambipolar diffusion calculations (e.g. Fiedler & Mouschovias 
1993). However, if turbulent motions in clouds are constrained to be compara- 
ble to the Alfven velocity, v a oc Bn~ - 5 , this result is also expected (Myers & 
Goodman 1988; Bertoldi & McKee 1992). 

The magnetic field direction, projected on the plane of the sky, can be mea- 
sured because spinning, aspherical grains tend to align their spin axes with the 
magnetic field direction (see Lazarian et al 1997 for a list of mechanisms). Then 
the dust grains absorb and emit preferentially in the plane perpendicular to the 
field. Consequently, background star light will be preferentially polarized along 
B± and thermal emission from the grains will be polarized perpendicular to B± 
(Hildebrand 1988). Goodman (1996) has shown that the grains that polarize 
background starlight do not trace the field very deeply into the cloud, but maps 
of polarized emission at far-infrared, submillimeter (Schleuning 1998) and mil- 
limeter (Akeson & Carlstrom 1997, Rao et al 1998) wavelengths are beginning 
to provide maps of field direction deep into clouds. Line emission may also be 
weakly polarized under some conditions (Goldreich & Kylafis 1981), providing 
a potential probe of B± with velocity information. After many attempts, this 
effect has been detected recently (Greaves et al 1999). 

The ionization fraction (x e ) is determined by chemical analysis and has been 
discussed by van Dishoeck and Blake (1998). Theoretically, x e should drop from 
about 10 -4 near the outer edge of the cloud to about 10 -8 in interiors shielded 
from ultraviolet radiation. Observational estimates of x e are converging on values 
around 10~ 8 to 1CT 7 in cores (de Boisanger et al 1996, Caselli et al 1998, Williams 
et al 1998, Bergin et al 1999). 

3.5 Observational and Analytical Tools 

Having discussed how different physical conditions are probed, I will end this 
section with a brief summary of the observational and analytical tools that are 
used. Clearly, most information on physical conditions comes from observations 
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of molecular lines. Most of these lie at millimeter or submillimeter wavelengths, 
and progress in this field has been driven by the development of large single-dish 
telescopes operating at submillimeter wavelengths and by arrays of antennas op- 
erating interferometrically at millimeter wavelengths (Sargent & Welch 1993). 
The submillimeter capability has allowed the study of high- J levels for excitation 
analysis and increased sensitivity to dust continuum emission, which rises with 
frequency (S v oc v 2 or faster). Studies of millimeter and submillimeter emission 
from dust have been greatly enhanced recently with the development of cameras 
on single dishes, both at millimeter wavelengths (Kreysa 1992) and at submillime- 
ter wavelengths, with SHARC (Hunter et al 1996) and SCUBA (Cunningham et 
al 1994). Examples of the maps that these cameras are producing are the color 
plates showing the 1.3 mm emission from the p Ophiuchi region (Motte et al 
1998, Figure 1) and the 850 pm and 450 fim emission from the ridge in Orion 
(L1641) (Johnstone & Bally 1999, Figure 2). 

Interferometric arrays, operating at millimeter wavelengths, have provided un- 
precedented angular resolution (now better than 1") maps of both molecular line 
and continuum emission. They are particularly critical for separating the contin- 
uum emission from a disk and the envelope and for studying deeply embedded 
binaries (e.g. Looney et al 1997, Figure 3). Complementary information has 
been provided in the infrared, with near-infrared star-counting (Lada et al 1994, 
1999), near- infrared and mid-infrared spectroscopy of rovibrational transitions 
(e.g. Mitchell et al 1990, Evans et al 1991, Carr et al 1995, van Dishoeck et al 
1998), and far- infrared continuum and spectral line studies (e.g. Haas et al 1995). 
Early results from the Infrared Space Observatory can be found in Yun & Liseau 
(1998). 

The analytical tools for molecular cloud studies have grown gradually in so- 
phistication. Early studies assumed LTE excitation, an approximation that is 
still used in some studies of CO isotopomers, but it is clearly invalid for other 
species. Studies of excitation require solution of the statistical equilibrium equa- 
tions (Goldsmith 1972). Goldreich h Kwan (1974) pointed out that photon 
trapping will increase the average T ex and provided a way of including its effects 
that was manageable with the limited computer resources of that time: the Large 
Velocity Gradient (LVG) approximation. Tied originally to their picture of col- 
lapsing clouds, this approximation allowed one to treat the radiative transport 
locally. Long after the overall collapse scenario had been discarded, the LVG 
method has remained in use, providing a quick way to include trapping, at least 
approximately. In parallel, more computationally intensive codes were developed 
for microturbulent clouds, in which photons emitted anywhere in the cloud could 
affect excitation anywhere else (e.g. Lucas 1974). 

The microturbulent and LVG assumptions are the two extremes, and real 
clouds probably lie between. For modest optical depths, the conclusions of the 
two methods differ by factors of about 3, comparable to uncertainties caused by 
uncertain geometry (White 1977, Snell 1981). These methods are still useful in 
some situations, but they are gradually being supplanted by more flexible radia- 
tive transport codes, using either the Monte Carlo technique (Bernes 1979, Choi 
et al 1995, Park & Hong 1998, Juvela 1997) or A-iteration (Dickel & Auer 1994, 
Yates et al 1997, Wiesemeyer 1999). Some of these codes allow variations in the 
velocity, density, and temperature fields, non-spherical geometries, clumps, etc. 
Of course, increased flexibility means more free parameters and the need for more 
extensive observations to constrain them. 
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Similar developments have occurred in the area of dust continuum emission. 
Since stellar photons are primarily at wavelengths where dust is quite opaque, a 
radiative transport code is needed to compute dust temperatures as a function of 
distance from a stellar heat source (Egan et al 1988, Manske & Henning 1998). 
For clouds without embedded stars or protostars, only the interstellar radiation 
field heats the dust; Tjj can get very low (5-10 K) in centers of opaque clouds 
(Leung 1975). Embedded sources heat clouds internally; in clouds opaque to 
the stellar radiation, it is absorbed close to the source and reradiated at longer 
wavelengths. Once the energy is carried primarily by photons at wavelengths 
where the dust is less opaque, the temperature distribution relaxes to the optically 
thin limit (Doty & Leung 1994): 



where L is the luminosity of the source and q = 2/((3 + 4), assuming oc v@ . 



Many low-mass stars actually form in regions of high-mass star formation, where 
clustered formation is the rule (Elmegreen 1985, Lada 1992, McCaughrean &; 
Stauffer 1994). The focus here is on regions where we can isolate the individual 
star-forming events and these are almost inevitably forming low-mass stars. 

4-1 Theoretical Issues 

The theory of isolated star formation has been developed in some detail. It relies 
on the existence of relatively isolated regions of enhanced density that can collapse 
toward a single center, though processes at smaller scales may cause binaries or 
multiples to form. One issue then is whether isolated regions suitable for forming 
individual, low-mass stars are clearly identifiable. The ability to separate these 
from the rest of the cloud underlies the distinction between sterile and fertile 
parts of clouds (§|]). 

Because of the enormous compression needed, gravitational collapse plays a key 
role in all star formation theories. In most cases, only a part of the cloud collapses, 
and theories differ on how this part is distinguished from the larger cloud. Is 
it brought to the verge of collapse by an impulsive event, like a shock wave 
(Elmegreen & Lada 1977) or a collision between clouds or clumps (Loren 1976), 
or is the process gradual? Among gradual processes, the decay of turbulence 
and ambipolar diffusion are leading contenders. If the decay of turbulence leaves 
the cloud in a subcritical state (M < Mb), then a relatively long period of 
ambipolar diffusion is needed before dynamical collapse can proceed (§|2|). If the 
cloud is supercritical (M > Mb), then the magnetic field alone cannot stop the 
collapse (e.g. Mestel 1985). If turbulence does not prevent it, a rapid collapse 
ensues, and fragmentation is likely. Shu et al (1987a, 1987b) suggested that the 
subcritical case describes isolated low mass star formation, while the supercritical 
case describes high-mass and clustered star formation. Recently, Nakano (1998) 
has argued that star formation in subcritical cores via ambipolar diffusion is 
implausible; instead he favors dissipation of turbulence as the controlling factor. 
For the present section, the questions are whether there is evidence that isolated, 
low-mass stars form in sub-critical regions and what is the status of turbulence 
in these cores. 
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Rotation could in principle support clouds against collapse, except along the 
rotation axis (Field 1978). Even if rotation does not prevent the collapse, it is 
likely to be amplified during collapse, leading at some point to rotation speeds 
able to affect the collapse. In particular, rotation is usually invoked to produce 
binaries or multiple systems on small scales. What do we know about rotation 
rates on large scales and how the rotation is amplified during collapse? Is there 
any correlation between rotation and the formation of binaries observationally? 
If rotation controls whether binaries form, can we understand why collapse leads 
to binary formation roughly half the time? It is clear that both magnetic flux and 
angular momentum must be redistributed during collapse to produce stars with 
reasonable fields and rotation rates, and these processes will affect the formation 
of binaries and protoplanetary disks. Some of these questions will be addressed 
in the next sections on globules and cores, while others will be discussed in the 
context of testing specific theories. 

4-2 Globules and Cores: Overall Properties 

Nearby small dark clouds, or globules, are natural places to look for isolated star 
formation (Bok & Reilly 1947). A catalog of 248 globules (Clemens & Barvainis 
1988) has provided the basis for many studies. Yun & Clemens (1990) found 
that 23% of the CB globules appear to contain embedded infrared sources, with 
spectral energy distributions typical of star forming regions (Yun 1993). About 
one-third of the globules with embedded sources have evidence of outflows (Yun 
& Clemens 1992; Henning & Launhardt 1998). Clearly, star formation does occur 
in isolated globules. 

Within the larger dark clouds, one can identify numerous regions of high opac- 
ity (e.g. Myers et al 1983), commonly called cores (Myers 1985). Surveys of 
such regions in low-excitation lines of NH3 (e.g. Benson & Myers 1989) led to 
the picture of an isolated core within a larger cloud, which then might pursue 
its course toward star formation in relative isolation from the rest of the cloud. 
Most intriguing was the fact that the NH3 linewidths in many of these cores in- 
dicated that the turbulence was subsonic (Myers 1983); in some cores, thermal 
broadening of NH3 lines even dominated over turbulent broadening (Myers & 
Benson 1983, Fuller & Myers 1993). Although later studies in other lines indi- 
cated a more complex dynamical situation (Zhou et al 1989, Butner et al 1995), 
the NH3 data provided observational support for theories describing the collapse 
of isothermal spheres (Shu 1977). The discovery of IRAS sources in half of these 
cores (Beichman et al 1986) indicated that they were indeed sites of star forma- 
tion. The observational and theoretical developments were synthesized into an 
influential paradigm for low mass star formation (Shu et al 1987a, § |4.5| ). 

Globules would appear to be an ideal sample for measuring sizes since the 
effects of the environment are minimized by their isolation, but distances are 
uncertain. Based on the angular size of the optical images and an assumed 
average distance of 600 pc (Clemens & Barvainis 1988), the mean size, (I) = 0.7 
pc. A subsample of these with distance estimates were mapped in molecular 
lines, yielding much smaller average sizes: (/) = 0.33 ± 0.15 pc for a sample of 
6 "typical" globules mapped in CS (Launhardt et al 1998); maps of the same 
globules in C 18 J = 2 — > 1 (Wang et al 1995) give sizes smaller by a factor of 
2.9 ± 1.6. A sample of 11 globules in the southern sky, mapped in NH3 (Bourke 
et al 1995) have (I) = 0.21 ± 0.08 pc. 
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Cores in nearby dark clouds have the advantage of having well-determined 
distances; the main issue is how clearly they stand out from the bulk of the 
molecular cloud. Gregersen (1998) found that some known cores are barely visible 
above the general cloud emission in the C 18 J = 1 — > line. The mean size of 
a sample of 16 cores mapped in NH3 is 0.15 pc, whereas CS J = 2 — > 1 gives 0.27 
pc, and C 18 J = 1 — > gives 0.36 pc (Myers et al 1991). These differences may 
reflect the effects of opacity, chemistry, and density structure. 

Globules are generally not spherical. By fitting the opaque cores with ellipses, 
Clemens k, Barvainis (1988) found a mean aspect ratio (a/b) of 2.0. Measure- 
ments of aspect ratio in tracers of reasonably dense gas toward globules give 
{a/b) ~ 1.5 — 2 (Wang et al 1995, Bourke et al 1995), and cores in larger clouds 
have (a/b) ~ 2 (Myers et al 1991). Myers et al and Ryden (1996) have argued 
that the underlying 3-D shapes were more likely to be prolate, with axial ratios 
around 2, than oblate, where axial ratios of 3-10 were needed. However, toroids 
may also be able to match the data because of their central density minima (Li 
& Shu 1996). 

The uncertainties in size are reflected directly into uncertainties in mass. For 
the sample of globules mapped by Bourke et al (1995) in NH3, (M^) = 4± 1 M , 
compared to 10 ±2 M for cores. Larger masses are obtained from other tracers. 
For the sample of globules studied by Launhardt et al (1998), (My) = 26 ± 12 
based on CS J = 2 -» 1 and 10 ± 6, based on C 18 J = 2 — > 1. Studies of the 
cores in larger clouds found similar ranges and differences among tracers (e.g. 
Fuller 1989, Zhou et al 1994a). A series of studies of the Taurus cloud complex 
has provided an unbiased survey of cores with known distance identified by C 18 
J = 1 — > maps. Starting from a large-scale map of 13 CO J = 1 — > (Mizuno 
et al 1995), Onishi et al (1996) covered 90% of the area with N > 3.5 x 10 21 
with a map of C 18 J = 1 — > with 0.1 pc resolution. They identified 40 cores 
with (I) = 0.46 pc, (a/b) = 1.8, and (Mjv) = 23 Mq. The sizes extended over a 
range of a factor of 6 and masses over a factor of 80. Comparing these cores to 
the distribution of T Tauri stars, infrared sources, and H 13 CO + emission, Onishi 
et al (1998) found that all cores with N > 8 x 10 21 cm~ 2 are associated with 
H 13 CO + emission and/or cold IRAS sources. In addition, the larger cores always 
contained multiple objects. They concluded that the core mass per star- forming 
event is relatively constant at 11 M . 

It is clear that characterizing globules and cores by a typical size and mass 
is an oversimplification. First, they come in a range of sizes that is probably 
just the low end of the general distribution of cloud sizes. Second, the size and 
mass depend strongly on the tracer and method used to measure them. If we 
ignore these caveats, it is probably fair to say that most of these regions have 
sizes measured in tracers of reasonably dense gas in the range of a few tenths 
of a pc and masses of less than 100 M , with more small, low mass cores than 
massive ones. The larger cores tend to be fragmented, so that the mass of gas 
with n > 10 4 cm -3 tends toward 10 M , within a factor of 2. Star formation 
may occur when the column density exceeds 8 x 10 21 cm~ 2 , corresponding in this 
case to n ~ 10 4 cm^ 3 . 

4-3 Globules and Cores: Internal Conditions 

If the sizes and masses of globules and cores are poorly defined, at least the 
temperatures seem well understood. Early CO observations (Dickman 1975) 
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showed that the darker globules are cold 10 K), as expected for regions with 

only cosmic ray heating. Similar results were found from NH3 (Martin & Barrett 
1978 Myers & Benson 1983; Bourke et al 1995). Clemens et al (1991) showed the 
distribution of CO temperatures for a large sample of globules; the main peak 
corresponded to Tk = 8.5 K, with a small tail to higher Tk- Determination of the 
dust temperature was more difficult, but Keene (1981) measured To = 13 — 16 K 
in B133, a starless core. More recently, Ward-Thompson et al (1998) used ISO 
data to measure To = 13 K in another starless core. These results are similar 
to predictions for cores heated by the interstellar radiation field (Leung 1975), 
though Td is expected to be lower in the deep interiors. 

A sample of starless globules with Ay ~ 1 — 2 mag produced no detec- 
tions of NH 3 (Kane et al 1994), but surveys of H 2 CO (Wang et al 1995) and 
CS (Launhardt et al 1998, Henning & Launhardt 1998) toward more opaque 
globules indicate that dense gas is present in some. The detection rate of CS 
J = 2 — > 1 emission was much higher in globules with infrared sources. Butner 
et al (1995) analyzed multiple transitions of DCO + in 18 low- mass cores, finding 
(log ra(cm -3 )) ~ 5, with a tendency to slightly higher values in the cores with 
infrared sources. Thus, gas still denser than the (logn(cm -3 )) ~ 4 gas traced by 
NH3 exists in these cores. 

Discussion of the kinematics of globules and cores has often focused on the 
relationship between linewidth (Av) and size (/ or R). This relation is much less 
clearly established for cores within a single cloud than is the relation for clouds 
as a whole (§ |3.3| ), and it may have a different origin (Myers 1985). Goodman 
et al (1998) have distinguished four types of linewidth-size relationships. Most 
studies have employed Goodman Type 1 relations (multi-tracer, multi-core), but 
it is difficult to distinguish different causes in such relations. Goodman Type 2 
relations (single-tracer, multi-core) relations within a single cloud would reveal 
whether the virial masses of cores are reliable. Interestingly, the most systematic 
study of cores in a single cloud found no correlation in 24 cores in Taurus mapped 
in C 18 J = 1 -> (7 = 0.0 ± 0.2), but the range of sizes (0.13 to 0.4 pc) may 
have been insufficient (Onishi et al 1996). 

To study the kinematics of individual cores, the most useful relations are Good- 
man Types 3 (multi-tracer, single-core) and 4 (single-tracer, single-core). In 
these, a central position is defined, either by an infrared source or a line peak. 
A Type 3 relationship using NH3, C 18 J = 1 — > and CS J = 2 — > 1 lines was 
explored by Fuller & Myers (1992), who found Av cx R? , with R the radius of 
the half-power contour. Caselli & Myers (1995) added 13 CO J = 1 — > data and 
constructed a Type 1 relation for 8 starless cores, after removing the thermal 
broadening. The mean 7 was 0.53 ± 0.07 with a correlation coefficient of 0.81. 
However, both these relationships depend strongly on the fact that the NH3 lines 
have small Av and R. Some other species (e.g. HC3N, Fuller &z Myers 1993) 
also have narrow lines and small sizes, but DCO + emission has much wider lines 
over a similar map size (Butner et al 1995), raising the possibility that chemical 
effects cause different molecules to trace different kinematic regimes within the 
same overall core. Goodman et al (1998) suggest that the DCO + is excited in 
a region outside the NH3 region and that the size of the DCO + region is un- 
derestimated. On the other hand, some chemical simulations indicate that NH3 
will deplete in dense cores while ions like DCO + will not (Rawlings et al 1992), 
suggesting the opposite solution. 

A way to avoid such effects is to use a Type 4 relation, searching for a correlation 
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between Avnt as spectra are averaged in rings of larger size, though line of sight 
confusion cannot be avoided, as one gives up the ability to tune the density 
sensitivity. This method has almost never been applied to tracers of dense gas, 
but Goodman et al (1998) use an indirect method to obtain a Type 4 relation 
for three clouds in OH, C 18 0, and NH3. The relations are very flat (7 = 0.1 to 
0.3), and 7 / has statistical significance only for OH, which traces the least 
dense gas. In particular, NH3 shows no significant Type 4 relation, having narrow 
lines on every scale (see also Barranco & Goodman 1998). Goodman et al (1998) 
interpret these results in terms of a "transition to coherence" at the scale of 0.1- 
0.2 pc from the center of a dense core. Inside that radius, the turbulence becomes 
subsonic and no longer decreases with size (Barranco & Goodman 1998). While 
this picture accords nicely with the idea that cores can be distinguished from the 
surroundings and treated as "units" in low-mass star formation, the discrepancy 
between values of Av measured in different tracers of the dense core (cf Butner 
et al 1995) indicate that caution is required in interpreting the NH3 data. 

Rotation can be detected in some low-mass cores, but the ratio of rotational to 
gravitational energy has a typical value of 0.02 on scales of 0.1 pc (Goodman et al 
1993). The inferred rotation axes are not correlated with the orientation of cloud 
elongation, again suggesting that rotation is not dynamically important on the 
scale of 0.1 pc. Ohashi et al (1997a) find evidence in several star-forming cores 
for a transition at r ~ 0.03 pc, inside of which the specific angular momentum 
appears to be constant at about 10 -3 km s _1 pc down to scales of about 200 AU. 

Knowledge of the magnetic field strength in cores and globules would be ex- 
tremely valuable in assessing whether they are subcritical or supercritical. Unfor- 
tunately, Zeeman measurements of these regions are extremely difficult because 
they must be done with emission, unless there is a chance alignment with a back- 
ground radio source. Crutcher et al (1993) detected Zeeman splitting in OH in 
only one of 12 positions in nearby clouds. Statistical analysis of the detection and 
the upper limits, including the effects of random orientation of the field, led to 
the conclusion that the data could not falsify the hypothesis that the clouds were 
subcritical. Another problem is that the OH emission probes relatively low den- 
sities (n ~ 10 3 cm -3 ). Attempts to use CN to probe denser gas produced upper 
limits that were less than expected for subcritical clouds, but the small sample 
size and other uncertainties again prevented a definitive conclusion (Crutcher et 
al 1996). Improved sensitivity and larger samples are crucial to progress in this 
field. At present, no clear examples of subcritical cores have been found (Crutcher 
1999a), but uncertainties are sufficient to allow this possibility. 

Onishi et al (1996) found that the major axis of the cores they identified in 
Taurus tended to be perpendicular to the optical polarization vectors and hence 
B±. Counterexamples are known in other regions (e.g. Vrba et al 1976) and 
the fact that optical polarization does not trace the dense portions of clouds 
(Goodman 1996) suggests that this result be treated cautiously. Further studies 
of B± using dust emission at long wavelengths are clearly needed. 

The median ionization fraction of 23 low-mass cores is 9 x 10~ 8 , with a range 
of logx e of —7.5 to —6.5, with typical uncertainties of a factor of 0.5 in the log 
(Williams et al 1998). Cores with stars do not differ significantly in x e from 
cores without stars, consistent with cosmic ray ionization. For a cloud with 
n = 10 4 cm -3 , the ambipolar diffusion timescale, tAD ~ 7 x 10 6 yr ~ 20i/j. If 
the cores are subcritical, they will evolve much more slowly than a free fall time. 
Recent comparisons of the line profiles of ionized and neutral species have been 
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interpreted as setting an upper limit on the ion-neutral drift velocity of 0.03 
km s _1 , consistent with that expected from ambipolar diffusion (Benson et al 
1998). 

To summarize the last two subsections, there is considerable evidence that 
distinct cores can be identified, both as isolated globules and as cores within 
larger clouds. While there is a substantial range of properties, scales of 0.1 pc in 
size and 10 M in mass seem common. The cores are cold (Tk ~ Tp ~ 10 K) 
and contain gas with n ~ 10 4 cm~ 3 , extending in some cases up to n ~ 10 5 cm -3 . 
While different molecules disagree about the magnitude of the effect, these cores 
seem to be regions of decreased turbulence compared to the surroundings. While 
no clear cases of subcritical cores have been found, the hypothesis that low- 
mass stars form in subcritical cores cannot be ruled out observationally. How 
these cores form is beyond the scope of this review, but again one can imagine 
two distinct scenarios: ambipolar diffusion brings an initially subcritical core 
to a supercritical state; or dissipation of turbulence plays a similar role in a 
core originally supported by turbulence (Myers & Lazarian 1998). In the latter 
picture, cores may build up from accretion of smaller diffuse elements (Kuiper et 
al 1996), perhaps the structures inferred by Falgarone et al (1998). 

4-4 Classification of Sources and Evolutionary Scenarios 

The IRAS survey provided spectral energy distributions over a wide wavelength 
range for many cores (e.g. Beichman et al 1986), leading to a classification scheme 
for infrared sources. In the original scheme (Lada & Wilking 1984, Lada 1987), 
the spectral index between 2 /xm and the longest observed wavelength was used 
to divide sources into three Classes, designated by roman numerals, with Class I 
indicating the most emission at long wavelengths. These classes rapidly became 
identified with stages in the emerging theoretical paradigm (Shu et al 1987a): 
Class I sources are believed to be undergoing infall with simultaneous bipolar 
outflow, Class II sources are typically visible T Tauri stars with disks and winds, 
and Class III sources have accreted or dissipated most of the material, leaving a 
pre-main-sequence star, possibly with planets (Adams et al 1987; Lada 1991). 

More recently, submillimeter continuum observations have revealed a large 
number of sources with emission peaking at still longer wavelengths. Some of 
these new sources also have infrared sources and powerful bipolar outflows, in- 
dicating that a central object has formed; these have been designated Class 
(Andre et al 1993). Andre & Montmerle (1994) argued that Class sources rep- 
resent the primary infall stage, in which there is still more circumstellar than 
stellar matter. Outflows appear to be most intense in the earliest stages, de- 
clining later (Bontemps et al 1996). Other cores with submillimeter emission 
have no IRAS sources and probably precede the formation of a central object. 
These were found among the "starless cores" of Benson &: Myers (1989), and 
Ward-Thompson et al (1994) referred to them as "pre-protostellar cores." The 
predestination implicit in this name has made it controversial, and I will use the 
less descriptive (and somewhat tongue-in-cheek) term, Class —1. There has also 
been some controversy over whether Class sources are really distinct from Class 
I sources or just more extreme versions. The case for Class sources as a distinct 
stage can be found in Andre et al (2000). 

While classification has an honored history in astronomy, serious tests of the- 
ory are facilitated by continuous variables. Myers & Ladd (1993) suggested that 
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we characterize the spectral energy distribution by the flux-weighted mean fre- 
quency, or more suggestively, by the temperature of a black body with the same 
mean frequency. The latter (Tj, i) was calculated by Chen et al (1995) for many 
sources and the following boundary lines in T\, \ were found to coincide with the 
traditional classes: < 70K for Class 0, 70 < < 650K for Class I, and 
650 < Tboi < 2800K for Class II. In a crude sense, T^i captures the "coolness" of 
the spectral energy distribution, which is related to how opaque the dust is, but 
it can be affected strongly by how much mid-infrared and near-infrared emission 
escapes, and thus by geometry. Other measures, such as the ratio of emission 
at a submillimeter wavelength to the bolometric luminosity {L smm / L^ i), may 
also be useful. One of the problems with all these measures is that the bulk of 
the energy for classes earlier than II emerges at far- infrared wavelengths, where 
resolution has been poor. Maps of submillimeter emission are showing that many 
near-infrared sources are displaced from the submillimeter peaks and may have 
been falsely identified as Class I sources (Motte et al 1998). Ultimately, higher 
spatial resolution in the far-infrared will be needed to sort out this confusion. 

4-5 Detailed Theories 

The recent focus on the formation of isolated, low-mass stars is at least partly 
due to the fact that it is more tractable theoretically than the formation of 
massive stars in clusters. Shu (1977) argued that collapse begins in a centrally 
condensed configuration and propagates outward at the sound speed (a); matter 
inside rj„j = at is infalling after a time t. Because the collapse in this model 
is self-similar, the structure can be specified at any time from a single solution. 
This situation is called inside-out collapse. In Shu's picture, the pre-collapse 
configuration is an isothermal sphere, with n(r) oc r~ p and p = 2. Calculations 
of core formation via ambipolar diffusion gradually approach a configuration with 
an envelope that is close to a power law, but with a core of ever-shrinking size 
and mass where p ~ (e.g. Mouschovias 1991). It is natural to identify the 
ambipolar diffusion stage with the Class —1 stage, and the inside-out collapse 
with Class to Class I. As collapse proceeds, the material inside rj„/ becomes 
less dense, with a power law approaching p = 1.5 in the inner regions, after a 
transition region where the density is not a power law. In addition, v(r) oc r~ ' 5 
at small r. 

There are many other solutions to the collapse problem (Hunter 1977, Foster 
and Chevalier 1993) ranging from the inside-out collapse to overall collapse (Lar- 
son 1969, Penston 1969). Henriksen et al (1997) have argued that collapse begins 
before the isothermal sphere is established; the inner {p = 0) core undergoes a 
rapid collapse that they identify with Class sources. They suggest that Class 
I sources represent the inside-out collapse phase, which appears only when the 
wave of infall has reached the p = 2 envelope. 

Either rotation or magnetic fields will break spherical symmetry. Terebey et 
al (1984) added slow rotation of the original cloud at an angular velocity of Q 
to the inside-out collapse picture, resulting in another characteristic radius, the 
point where the rotation speed equals the infall speed. A rotationally supported 
disk should form somewhere inside this centrifugal radius (Shu et al 1987a) 
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where M is the mass already in the star and disk. Since disks are implicated 
in all models of the ultimate source of the outflows, the formation of the disk 
may also signal the start of the outflow. Once material close to the rotation axis 
has accreted (or been blown out) all further accretion onto the star should occur 
through the disk. 

Core formation in a magnetic field should produce a flattened structure (e.g. 
Fiedler & Mouschovias 1993) and Li and Shu (1996) argue that the equilibrium 
structure equivalent to the isothermal sphere is the isothermal toroid. Useful 
insights into the collapse of a magnetized cloud have resulted from calculations 
in spherical geometry (Safier et al 1997, Li 1998) or for thin disks (Ciolek & 
Konigl 1998). Some calculations of collapse in two dimensions have been done 
(e.g. Fiedler & Mouschovias 1993). A magnetically-channeled, flattened structure 
may appear; this has been called a pseudo-disk (Galli & Shu 1993a, 1993b) to 
distinguish it from the rotationally supported disk. In this picture, material would 
flow into a pseudo-disk at a scale of ~ 1000 AU before becoming rotationally 
supported on the scale of r c . The breaking of spherical symmetry on scales of 
1000 AU may explain some of the larger structures seen in some regions (see 
Mundy et al 2000 for a review). 

Ultimately, theory should be able to predict the conditions that lead to binary 
formation, but this is not presently possible. Steps toward this goal can be seen 
in numerical calculations of collapse with rotation (e.g. Bonnell & Bastien 1993, 
Truclove et al 1998, Boss 1998). Rotation and magnetic fields have been combined 
in a series of calculations by Basu & Mouschovias (1995 and references therein); 
Basu (1998) has considered the effects of magnetic fields on the formation of 
rotating disks. 

Theoretical models make predictions that are testable, at least in principle. In 
the simplest picture of the inside-out collapse of the rotationless, non-magnetic 
isothermal sphere, the theory predicts all the density and velocity structure with 
only the sound speed as a free parameter. Rotation adds O and magnetic fields 
add a reference field or equivalent as additional parameters. Departure from 
spherical symmetry adds the additional observational parameter of viewing angle. 
Consequently, observational tests have focused primarily on testing the simplest 
models. 

4-6 Tests of Evolutionary Hypotheses and Theory 

Both the empirical evolutionary sequence based on the class system and the 
detailed theoretical predictions can be tested by detailed observations. One can 
compute the expected changes in the continuum emission as a function of time 
for a particular model. Examples include plots of L versus Ay (Adams 1990), 
L versus T^ (Myers et al 1998), and L mm versus L (Saraceno et al 1996). For 
example, as time goes on, Ay decreases, T^i increases, and L reaches a peak 
and declines. At present the models are somewhat idealized and simplified, and 
different dust opacities need to be considered. Comparison of the number of 
objects observed in various parts of these diagrams with those expected from 
lifetime considerations can provide an overall check on evolutionary scenarios 
and provide age estimates for objects in different classes (e.g. Andre et al 2000). 
One can also test the models against observations of particular objects. However, 
models of source structure constrained only by the spectral energy distribution 
are not unique (Butner et al 1991, Men'shchikov & Henning 1997). Observational 
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determination of n(r) and v(r) can apply more stringent tests. 

Maps of continuum emission from dust can trace the column density as a func- 
tion of radius quite effectively, if the temperature distribution is known. With an 
assumption about geometry, this information can be related to n(r). For the Class 
— 1 sources, with no central object, the core should be isothermal (To ~ 5 — 10 
K) or warmer on the outside if exposed to the interstellar radiation field (Leung 
1975, Spencer & Leung 1978). New results from ISO will put tighter limits on 
possible internal energy sources (e.g. Ward-Thompson et al 1998, Andre et al 
2000). Based on small maps of submillimeter emission, Ward-Thompson et al 
(1994) found that Class —1 sources were not characterized by single power laws 
in column density; they fit the distributions with broken power laws, indicating 
a shallower distribution closer to the center. This distribution appears consistent 
with the interpretation that these are cores still forming by ambipolar diffusion. 
Maps at 1.3 mm (Andre et al 1996, Ward-Thompson et al 1999) confirm the 
early results: the column density is quite constant in an inner region (r < 3000 
to 4000 AU), with M ~ 0.7 M Q . SCUBA maps of these sources are just becom- 
ing available, but they suggest similar conclusions (D Ward-Thompson, personal 
communication, Shirley et al 1998). In addition, some of the cores are filamentary 
and fragmented (Ward-Thompson et al 1999). Citing these observed properties, 
together with a statistical argument regarding lifetimes, Ward-Thompson et al 
(1999) now argue that ambipolar diffusion in subcritical cores does not match 
the data (see also Andre et al 2000). 

For cores with central sources (Class > —1), Tn{r) will decline with radius 
from the source. If the emission is in the Rayleigh-Jeans limit, I u oc Te>k(v)N . 
If k(v) is not a function of r, I v {9) oc f T£>(r)n(r)dl, where the integration is 
performed along the line of sight (dl). If one avoids the central beam and any 
outer cut-off, assumes the optically thin expression for Tr,{r) oc r~ q , and fits a 
power law to I v {6) oc 6~ m , then p = m + 1 — q (Adams 1991, Ladd et al 1991). 
With a proper calculation of To(r) and convolution with the beam, one can 
include all the information (Adams 1991). This technique has been applied in 
the far-infrared (Butner et al 1991), but it is most useful at longer wavelengths. 
Ladd et al (1991) found p = 1.7±0.3 in two cores, assuming q = 0.4, as expected 
for (3 = 1 (equation 14). Results so far are roughly consistent with theoretical 
models, but results in this area from the new submillimeter cameras will explode 
about the time this review goes to press. 

Another issue arises for cores with central objects. If a circumstellar disk 
contributes significantly to the emission in the central beam, it will increase the 
fitted value of m. Disks contribute more importantly at longer wavelengths, so 
the disk contribution is more important at millimeter wavelengths (e.g. Chandler 
et al 1995). Luckily, interferometers are available at those wavelengths, and 
observations with a wide range of antenna spacings can separate the contributions 
of a disk and envelope. Application of this technique by Looney et al (1997) to 
L1551 IRS5, a Class I source, reveals very complex structure (Figure 3): binary 
circumstellar disks (cf Rodriguez et al 1998), a circumbinary structure (perhaps 
a pseudo-disk), and an envelope with a density distribution consistent with a 
power law (p = 1.5 to 2). This technique promises to be very fruitful in tracing 
the flow of matter from envelope to disk. Early results indicated that disks 
are more prominent in more evolved (Class II) systems (Ohashi et al 1996), 
but compact structures are detectable in some younger systems (Terebey et al 
1993). Higher resolution observations and careful analysis will be needed to 
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distinguish envelopes, pseudo-disks, and Keplerian disks (see Mundy et al 2000 
for a review). At the moment, one can only say that disks in the Class stage 
are not significantly more massive than disks in later classes (Mundy et al 2000). 
Meanwhile, the interferometric data confirm the tendency of envelope mass to 
decrease with class number inferred from single-dish data (Mundy et al 2000). 



L1551 IRS5 A = 2.7 mm Visibilities 
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Figure 3: Observed visibilities of L1551 IRS 5 at 2.7 mm, binned in annuli (open 
squares with error bars), are plotted versus the projected baseline in units of 10 3 
times the wavelength. Different models are shown in each panel by the small 
boxes. Panel a has a model with only a Gaussian source of radius 80 AU; panel b 
has a model with two point sources constrained to match the map; panel c adds 
a truncated power law (p = 1.5, q = 0.5) envelope to the two point sources; panel 
d adds to the previous components a circumbinary structure, represented by a 
Gaussian. In the final optimization (panel d), the envelope has a mass of 0.28 
M and an outer radius of 1100 AU, the circumbinary structure has a mass of 
0.04 M Q , and the circumstellar disk masses are 0.024 and 0.009 M (Looney et 
al 1997). 

Similar techniques have been used for maps of molecular line emission. For 
example, 13 CO emission has been used to trace column density in the outer 
regions of dark clouds. With an assumption of spherical symmetry, the results 
favor p ~ 2 in most clouds (Snell 1981, Arquilla k Goldsmith 1985). The 13 CO 
lines become optically thick in the inner regions; studies with higher spatial 
resolution in rarer isotopomers, like C 18 or C 17 0, tend to show somewhat more 
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shallow density distributions than expected by the standard model (Zhou et al 
1994b, Wang et al 1995). Depletion in the dense, cold cores may still confuse 
matters (e.g. Kuiper et al 1996, §3.1). Addressing the question of evolution, 
Ladd et al (1998) used two transitions of C 18 and C 17 to show that N toward 
the central source declines with T^oi, with a power between 0.4 and 1.0. To 
reproduce the inferred rapid decrease in mass with time, they suggest higher 
early mass loss than predicted by the standard model. 

By observing a series of lines of different critical density, modeling those lines 
with a particular cloud model and appropriate radiative transport, and predicting 
the emission into the beams used for the observations, one can constrain the 
run of density more directly. Studies using two transitions of H2CO have again 
supported p = 2 ± 0.5 on relatively large scales (Loren et al 1983, Fulkerson &; 
Clark 1984). When interferometery of H2CO was used to improve the resolution 
on one core, p appeared to decrease at small r (Zhou et al 1990), in agreement 
with the model of Shu (1977). Much of the recent work on this topic has involved 
testing of detailed collapse models, including velocity fields and the complete 
density law, rather than a single power law, as described in the next section. 



4-7 Collapse 

The calculation of line profiles as a function of time (Zhou 1992) for the collapse 
models of Shu (1977) and Larson (1969) and Penston (1969), along with claims of 
collapse in a low-mass star forming region (Walker et al 1986), reinvigorated the 
study of protostellar collapse. Collapsing clouds will depart from the linewidth- 
size relation (§[D|), having systematically larger linewidths for a given size (Zhou 
1992). Other simulations of line profiles range from a simple two-layer model 
(Myers et al 1996) to detailed calculations of radiative transport (Choi et al 
1995, Walker et al 1994, Wiesemeyer 1997, 1999). 

Zhou et al (1993) showed that several lines of CS and H2CO observed towards 
B335, a globule with a Class source, could be fitted very well using the exact n(r) 
and v(r) of the inside-out collapse model. Using a more self-consistent radiative 
transport code, Choi et al (1995) found slightly different best-fit parameters. 
Using a sound speed determined from lines away from the collapse region, the 
only free parameters were the time since collapse began and the abundance of each 
molecule. With several lines of each molecule, the problem is quite constrained 
(Figure 4). This work was important in gaining acceptance for the idea that 
collapse had finally been seen. 

Examination of the line profiles in Figure 4 reveals that most are strongly self- 
absorbed. Recall that the overall collapse idea of Goldreich and Kwan (1974) 
was designed to avoid self-absorbed profiles. The difference is that Goldreich and 
Kwan assumed that v(r) cx r, so that every velocity corresponded to a single 
point along the line of sight. In contrast, the inside-out collapse model predicts 
v(r) oc r~ 5 inside a static envelope. If the line has substantial opacity in the 
static envelope, it will produce a narrow self-absorption at the velocity centroid 
of the core (Figure 5). The other striking feature of the spectra in Figure 4 is 
that the blue-shifted peak is stronger than the red-shifted peak. This "blue" 
profile occurs because the v(r) oc r™ 0,5 velocity field has two points along any 
line of sight with the same Doppler shift (Figure 6). For a centrally peaked 
temperature and density distribution, lines with high critical densities will have 
higher T ex at the point closer to the center. If the line has sufficient opacity at the 
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Figure 4: Line profiles of H2CO and CS emission (solid histograms) toward B335 
and the best-fitting model (dashed lines). The model line profiles were calculated 
with a Monte Carlo code, including non-LTE excitation and trapping, with an 
input density and velocity field taken directly from the collapse model of Shu 
(1977) and a temperature field calculated with a separate dust radiation transport 
code. The best-fitting model has an infall radius of 0.03 pc (Choi et al 1995). 

relevant point in the cloud, the high T ex point in the red peak will be obscured 
by the lower T ex one, making the red peak weaker than the blue peak (Figure 
6). Thus a collapsing cloud with a velocity and density gradient similar to those 
in the inside-out collapse model will produce blue profiles in lines with suitable 
excitation and opacity properties. A double-peaked profile with a stronger blue 
peak or a blue-skewed profile relative to an optically thin line then becomes a 
signature of collapse. These features were discussed by Zhou & Evans (1994) and, 
in a more limited context, by Snell & Loren (1977) and Leung & Brown (1977). 

Of course, the collapse interpretation of a blue profile is not unique. Such 
profiles can be produced in a variety of ways. To be a plausible candidate for 
collapse, a core must also show these features: an optically thin line must peak 
between the two peaks of the opaque line; the strength and skewness should 
peak on the central source; and the two peaks should not be caused by clumps 
in an outflow. The optically thin line is particularly crucial, since two cloud 
components, for example, colliding fragments, could produce the double-peaked 
blue profile, but they would also produce a double-peaked profile in the optically 
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Figure 5: The origin of various parts of the line profile for a cloud undergoing 
inside-out collapse. The static envelope outside ri n f produces the central self- 
absorption dip, the blue peak comes from the back of the cloud, and the red peak 
from the front of the cloud. The faster collapse near the center produces line 
wings, but these are usually confused by outflow wings. 

thin line. 

Rotation, combined with self-absorption, can create a line profile like that of 
collapse (Menten et al 1987, Adelson & Leung 1988), but toward the center of 
rotation, the line would be symmetric (Zhou 1995). Rotating collapse can cause 
the line profiles to shift from blue to red-skewed on either side of the rotation 
axis, with the sign of the effect depending on how the rotation varies with radius 
(Zhou 1995). Maps of the line centroid can be used to separate rotation from 
collapse (Adelson & Leung 1988, Walker et al 1994). 

To turn a collapse candidate into a believable case of collapse, one has to map 
the line profiles, account for the effects of outflows, model rotation if present, 
and show that a collapse model fits the line profiles. To date this has been 
done only for a few sources: B335 (Zhou et al 1993, Choi et al 1995), L1527 
(Myers et al 1995, Zhou et al 1996, Gregersen et al 1997), and IRAS 16293 (Zhou 
1995, Narayanan et al 1998). Of this group, only IRAS 16293, rotating about 
20 times faster than B335, is known to be a binary (Wootten 1989, Mundy et al 
1992), supporting the idea that faster rotation is more likely to produce a binary. 
Mathieu (1994) reviews binarity in the pre-main-sequence stage, and Mundy et 
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al (2000) discuss recent evidence on the earlier stages. 

Interferometric observations have also revealed infall motions and rotational 
motions on scales of ~ 1000 AU in several sources (e.g. Ohashi et al 1997b; 
Momose et al 1998). Such studies can reveal how matter makes the transition 
from infall to a rotating disk. Inevitably, irregularities in the density and velocity 
fields will confuse matters in real sources, and these may be more noticeable with 
interferometers. Outflows are particularly troublesome (Hogerheijde et al 1998). 
Extreme blue/red ratios are seen in interferometric observations of HCO + and 
HCN J = 1 — > lines, which are difficult to reproduce with standard models 
(Choi et al 1999). Even in B335, the best case for collapse, Velusamy et al (1995) 
found evidence for clumpy structure within the overall gradients. In addition, 
very high resolution observations of CS J = 5 — > 4 emission toward B335 are not 
consistent with predicted line profiles very close to the forming star (Wilner et al 
1999); either CS is highly depleted in the infalling gas, or the velocity or density 
fields depart from the model. 

If, for the sake of argument, we accept a blue profile as a collapsing core, can 
we see any evolutionary trends? A number of surveys for blue profiles have been 
undertaken recently. Gregersen et al (1997) found 9 blue profiles in 23 Class 
sources, using the J = 3 — > 2 and J = 4 — > 3 lines of HCO + and H 13 CO + . After 
consideration of possible confusion by outflows, etc., they identified 6 sources as 
good candidates for collapse. Mardones et al (1997) extended the search to Class 
I sources with T^i < 200 K, using CS and H2CO lines. They introduced the line 
asymmetry as a collapse indicator: 

5V = (V thlck - V thm )/AV thm , (16) 

where Vthick is the velocity of the peak of the opaque line, Vthin is the velocity 
of the peak of the optically thin line, and AVthin is the linewidth of the thin 
line. They confirmed many of the collapse candidates found by Gregersen and 
identified 6 more, but they found very few collapse candidates among the Class 
I sources. The difference could be caused by using different tracers, since the CS 
and H2CO lines are less opaque than the HCO + lines. To remove this uncer- 
tainty, Gregersen (1998) surveyed the Class I sources of Mardones et al (1997) in 
HCO + . Using SV as the measure, the fraction of blue profiles did not decrease 
substantially between Class and Class I. Most of these line profiles need further 
observations before they become bona fide candidates. 

When does collapse begin? Surveys of the Class —1 cores might reveal very 
early stages of collapse. In the inside-out collapse picture, blue profiles should 
appear, if at all, only toward the center. In fact, blue profiles have been seen in a 
substantial number of these cores (Gregersen 1998, Lee et al 1999) and maps of 
one, L1544, show that the blue profiles are very extended spatially (Tafalla et al 
1998). Clearly, extended blue profiles in Class —1 cores do not fit expectations 
for early collapse, and Tafalla et al argue that the velocities are higher than 
expected for ambipolar diffusion. If the regions producing the blue and red peaks 
are indeed approaching one another, they are forming a core more rapidly than in 
standard models, suggesting that some new ideas may be necessary (e.g. Nakano 
1998, Myers & Lazarian 1998). For a current review of this field, see Myers et al 
(2000). 
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Figure 6: A schematic explanation of why line profiles of optically thick, high- 
excitation lines are skewed to the blue in a collapsing cloud. The ovals are loci of 
constant line- of- sight velocity, for v(r) oc r~ ' 5 . Each line of sight intersects these 
loci at two points. The point closer to the center will have a higher T ex , especially 
in lines that are hard to excite, so that T ex (i?2) > T ex (Ri) and T ex (I?2) > T ex (B\). 
If the line is sufficiently opaque, the point R\ will obscure the brighter R2, but 
B2 lies in front of B\. The result is a profile with the blue peak stronger than 
the red peak (Zhou & Evans 1994). 

4-8 Summary of Isolated Star Formation 

Distinct cores can be identified in tracers of dense (n > 10 4 cm -3 ) gas; these 
cores are frequently associated with star formation. There is no clear evidence 
that they are magnetically subcritical, and some kinematic evidence suggests that 
the decay of turbulence, rather than ambipolar diffusion, is the critical feature. 
An empirical evolutionary sequence, based on the spectral appearance of dust 
emission, and detailed theoretical models are now being tested by observations. 
The spatial distribution of dust emission is providing important tests by probing 
n(r). Predictions of the evolution of spectral lines during collapse are available 
for the simplest theory and observations of some sources match predictions of 
theory quite well. Evidence of collapse is now strong in a few cases and surveys 
for distinctive line profiles have revealed many more possible cases. 
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5 CLUSTERED STAR FORMATION AND MASSIVE STARS 

In this section, we will address the issues of clustered formation, regardless of 
mass, and high-mass star formation, which seems to occur exclusively in clusters. 
In this review, the term cluster refers to a group of forming stars, whether or 
not the eventual outcome is a bound cluster. Since high mass stars are rare, 
the nearest examples are more distant than is the case for low mass star for- 
mation. Together with the fact that they form in clusters, the greater distance 
makes it difficult to isolate individual events of star formation. On the other 
side of the balance, massive stars are more easily detectable at large distances 
because luminosity is such a strong function of mass. Heating of the surroundings 
makes them strong emitters in both dust and many spectral lines; the spectra 
of regions forming massive stars are often very rich, reflecting both a high exci- 
tation state and, in many cases, enhanced abundances, as complex chemistry is 
driven by elevated temperatures (van Dishoeck & Blake 1998). These features led 
early studies, when sensitivity was poor, to concentrate on regions forming mas- 
sive stars. However, most of these advantages arise because the star is strongly 
influencing its surroundings; if we wish to know the preconditions, we will be 
misled. This problem is aggravated by the fast evolution of massive stars to the 
main sequence. Reviews of the topic of clustered star formation can be found in 
Elmegreen (1985), Lada & Lada (1991), and Lada (1999). Reviews focusing on 
the formation of massive stars include Churchwell (1993, 1999), Walmsley (1995), 
Stahler et al (2000), and Kurtz et al (2000). 

5. 1 Theoretical Issues 

Some of the primary theoretical issues regarding the formation of massive stars 
have been reviewed by Stahler et al (2000). First, what is the relevant dividing 
line between low-mass and high-mass stars? For the star formation problem, the 
question is how far in mass the scenario for low-mass stars can be extended. The- 
oretically, the limit is probably about 10 M 6 , where stars reach the main sequence 
before the surrounding envelope is dissipated. Observations of the physical con- 
ditions in regions forming intermediate mass stars (Herbig Ae/Be stars and their 
more embedded precursors) can reveal whether modifications are needed at even 
lower masses. Since accretion through disks plays a crucial role in the standard 
model, it is important to know the frequency and properties of disks around more 
massive stars. 

Stars as massive as 100 M Q seem to exist (Kudritzki et al 1992), but radia- 
tion pressure from the rapidly evolving stellar core should stop accretion before 
such masses can be built (e.g. Wolfire & Cassinelli 1987). In addition, mas- 
sive stars produce very strong outflows (Shepherd & Churchwell 1996, Bachiller 
1996). Since standard accretion theory cannot produce rates of mass accretion 
high enough to overwhelm these dispersive effects, some new effects must become 
important. 

Related questions concern the formation of clusters. To what extent can the 
ideas of isolated core collapse be applied if there are competing centers of collapse 
nearby? If the collapse to form massive stars is supercritical, a whole region may 
collapse and fragment to form many stars. 

The fact that the most massive stars are found near the centers of forming 
clusters has led to the suggestion that massive stars are built by collisional coa- 
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lescence of stars or protostars (Bonnell et al 1997, 1998). This scheme requires 
high densities (n* > 10 4 stars pc -3 ); for a mean stellar mass of 1 M , this corre- 
sponds to n > 2 x 10 5 cm" 3 . Even higher stellar densities are seen in the core of 
the Orion Nebula Cluster (Hillenbrand & Hartmann 1998). 

The special problems of making the most massive stars are a subset of the larger 
question of explaining the mass distribution of all stars. There may be variations 
in the IMF between clusters (Scalo 1998), which can test theories. Hipparcos 
observations of nearby OB associations have extended the membership to lower 
mass stars (de Zeeuw et al 1999), suggesting total masses of a few 10 3 M . 

The main questions these issues raise for observations are whether the mass 
and density of cores are sufficient for forming clusters and massive stars, and 
whether the mass distribution of clumps in a star forming region can be related 
to the mass distribution of stars ultimately formed. 



5.2 Overall Cloud and Core Properties 

What do we know about the general properties of the galactic clouds? The 
broadest picture is provided by surveys of CO and 13 CO. Surveys of significant 
areas of the Galaxy indicate that the power-law distribution in mass seen for 



small clouds ( dN{M) oc M~ a dM, § pTT|) continues up to a cutoff at M ~ 6 X 10 6 



M (Williams & McKee 1997). Studies by Casoli et al (1984), Solomon & Rivolo 
(1989), and Brand & Wouterloot (1995) find 1.4 < a < 1.8 over both inner and 
outer Galaxy. Extinction surveys find flatter slopes (Scalo 1985), but Scalo &; 
Lazarian (1996) suggest that cloud overlap affects the extinction surveys. The 
fact that a < 2 implies that most of the mass is in the largest structures, although 
there are issues of how to separate clouds at the largest scales. Since the star 
formation rate per unit mass, measured by CO, appears not to depend on cloud 
mass (Mead et al 1990, Evans 1991), the mass distribution supports the idea that 
most stars form in massive clouds (Elmegreen 1985). However, the enormous 
spread (> 10 2 ) in star formation rate per unit mass at any given mass, together 
with the fact that most of the molecular gas is sterile, suggests that comparisons 
to overall cloud masses, measured by CO, are not particularly relevant. 

Surveys in molecular lines indicative of denser gas have generally been biased 
towards signposts of star formation. Exceptions are the CS J = 2 — > 1 survey of 
L1630 (Lada et al 1991) and the CS J = 1 -> and J = 2 -> 1 surveys of L1641 
(Tatematsu et al 1993, 1998). These two clouds, also called Orion B and Orion 
A, are adjacent. In both cases, the CS J = 2 — > 1 maps showed more contrast 
than maps of 13 CO (Bally et al 1987), with CS J = 1 — > being somewhat in- 
termediate. Maps of higher- J transitions have been less complete, but show still 
less area covered by emission. The CS surveys detected less than 20% of the total 
mass in both clouds (Lada et al 1991, Tatematsu, personal communication). The 
J = 2 — > 1 emission in L1641 is somewhat smoother than the J = 2 — > 1 emis- 
sion from L1630 (Tatematsu et al 1998), and this difference may be reflected in 
the distribution of star formation. Star formation in L1641 appears to include a 
distributed component (Strom et al 1993); in contrast, star formation in L1630 is 
tightly concentrated in clusters associated with massive cores of dense gas (Lada 
1992, Li et al 1997). 

Excitation analysis of CS lines with higher critical density in L1630 shows 
that the star-forming regions all contain gas with n > 10 5 cm" 3 (Lada et al 
1997). These results suggest that surveys in lines of high n c (jk) are relevant for 
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characterizing star formation regions. Of the possible tracers, CS and NH3 have 
been most widely surveyed. In comparison to cores in Taurus, where only low- 
mass stars are forming, cores in the Orion clouds tend to be more massive and to 
have larger linewidths when observed with the same tracer (CS: Tatematsu et al 
1993, NH3: Harju et al 1993). The differences are factors of 2-4 for the majority 
of the cores, but larger for the cores near the Orion Nebula and those forming 
clusters in L1630 (Lada et al 1991). 

CS transitions have been surveyed toward Ultra-Compact (UC) HII regions, 
H2O masers, or luminous IRAS sources (see Kurtz et al 2000). Since the IRAS 
survey became available, most samples are drawn from the IRAS catalog with 
various color selection criteria applied. The most complete survey (Bronfman et 
al 1996) was toward IRAS sources with colors characteristic of UC HII regions 
(Wood and Churchwell 1989) over the entire Galactic plane. Bronfman et al 
found CS J = 2 — > 1 emission (see Table 1 for density sensitivity) in 59% of 1427 
IRAS sources, and the undetected sources were either weak in the far-infrared 
or had peculiar colors. Searches toward H2O masers have used the catalogs of 
Braz and Epchtein (1983) and Cesaroni et al (1988). Surveys of CS J = 2 — > 1 
(Zinchenko et al 1995, Juvela 1996) toward southern H2O masers found detection 
rates close to 100%, suggesting that dense, thermally excited gas surrounds the 
compact, ultra-dense regions needed to produce H2O masers. The detection rate 
drops in higher J transitions of CS (Plume et al 1992, 1997) but is still 58% 
in the CS J = 7 — > 6 line, which probes higher densities (Table 1). An LVG, 
multitransition study of CS lines found (logn(cm~ 3 )) = 5.9 for 71 sources and a 
similar result from a smaller sample using C 34 S data (Plume et al 1997). Densities 
derived assuming LVG fall between the average and maximum densities in clumpy 
models with a range of densities (Juvela 1997, 1998). 

Maps of the cores provide size and mass information. Based on the sizes and 
masses of 28 cores mapped in the CS J = 2 — > 1 line (Juvela 1996), one can 
compute a mean size, (I) = 1.2 ± 0.5 pc, mean virial mass, {My) ~ 5500 M , 
and (Afjv) ~ 4900 Mq. While the two mass estimates agree on average, there 
can be large differences in individual cases. Remarkably, cloud structure does not 
introduce a big uncertainty into the cloud masses: using a clumpy cloud model, 
(see § l5.4j ), Juvela (1998) found that M^j increased by a factor of 2 on average 
compared to homogeneous models and agreed with My to within a factor of 2. 
Plume et al (1997) obtained similar results from strip maps of CS J = 5 — > 4: 
(I) = 1.0 ± 0.7 pc (average over 25 cores); (My) = 3800 M (16 cores). As usual 
mean values must be regarded with caution; there is a size distribution. As cores 
with weaker emission are mapped, the mean size decreases; an average of 30 cores 
with full maps of J = 5 — > 4 emission gives (I) = 0.74 ± 0.56 pc, with a range of 
0.2 to 2.8 pc (Y Shirley, unpublished results). 

Churchwell et al (1992) surveyed 11 UC HII regions for CS J = 2 -> 1 and 
J = 5 — > 4 emission, leading to estimates of n > 10 5 cm~ 3 . Cesaroni et al (1991) 
surveyed 8 UC HII in three transitions of CS and C 34 S and estimated typical sizes 
of 0.4 pc, masses of 2000 Mq, and densities of 10 6 cm . More extensive surveys 
have been made in NH 3 ; Churchwell et al (1990) found NH 3 (J,K) = (1, 1) and 
(2,2) emission from 70% of a sample of 84 UC HII regions and IRAS sources 
with similar colors. They derived T%, finding a peak in the distribution around 
20 K, but a significant tail to higher values. Further studies in the (4,4) and 
(5,5) lines toward 16 UC HII regions with strong (2,2) emission (Cesaroni et al 
1992) detected a high fraction. Estimates for Tr ranged from 64 to 136 K and 
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sizes of 0.5 pc. Two sources indicated much higher densities and NH3 abundances. 
Follow-up studies with the VLA (Cesaroni et al 1994, 1998) revealed small (~ 0.1 
pc), hot (Tk ~ 50 — 200 K), dense (n = 10 7 cm -3 ) regions with enhanced NH3 
abundances. These hot cores (discussed below) are slightly displaced from the 
UC HII, but coincide with H2O masers. 

Magnetic fields strengths have been measured with the Zeeman effect toward 
about 10 regions of massive star formation. The fields are substantially stronger 
than the fields seen in isolated, low-mass cores, but the masses are also much 
higher. In most cases the mass to flux ratio is comparable to the critical ratio, 
once geometrical effects are considered (Crutcher 1999b). Given the uncertainties 
and sample size, it is too early to decide if regions forming massive stars are more 
likely to be supercritical than regions forming only low-mass stars. 

The ionization fraction in the somewhat more massive Orion cores appears to 
be very similar to that in low-mass cores: —6.9 < logx e < —7.3 (Bergin et al 
1999). The most massive cores in their sample have x e < 10 -8 , as do some of the 
massive cores studied by de Boisanger et al (1996). Expressed in terms of column 
density, the decline in x e appears around N ~ 3 x 10 22 cm -2 . At x e = 10~ 8 , 
tAD = 7 x 10 5 yr, about d.ltAD in isolated, low-mass cores. Even if the massive 
cores are subcritical, their evolution should be faster than that of low-mass cores. 

To summarize, the existing surveys show ample evidence that massive star 
formation usually takes place in massive (M > 10 3 M Q ), dense (n ~ 10 6 cm~ 3 ) 
cores, consistent with the requirements inferred from the study of young clusters 
and associations, and with conditions needed to form the most massive stars 
by mergers. Cores with measured B z seem to be near the boundary between 
subcritical and supercritical cores. 

5.3 Evolutionary Scenarios and Detailed Theories 

To what extent can an evolutionary scenario analogous to the class system be 
constructed for massive star formation? Explicit attempts to fit massive cores 
into the class system have relied on surveys of IRAS sources. Candidates for 
massive Class objects, with L > 10 3 L Q , have been found (Wilner et al 1995, 
Molinari et al 1998). One difficulty with using the shape of the spectral energy 
distribution for massive star formation is that dense, obscuring material usually 
surrounds objects even after they have formed, and a single star may be quite 
evolved but still have enough dust in the vicinity to have the same spectral 
energy distribution as a much younger object. The basic problem is the difficulty 
in isolating single objects. Also, the role of disks in massive regions is less clear, 
and they are unlikely to dominate the spectrum, as they do in low mass Class 
II sources. Other markers, such as the detection of radio continuum emission, 
must be used as age indicators. Hot cores provide obvious candidiates to be 
precursors of UC HII regions, but some have embedded UC HII regions and may 
be transitional (Kurtz et al 2000). The chemical state of massive cores may also 
provide an evolutionary sequence; Helmich et al (1994) suggested an evolutionary 
ordering of three sources in the W3 region based on their molecular spectra and 
chemical models (see van Dishoeck &: Blake 1998). 

While theories for clustered and massive star formation are much less devel- 
oped, some steps have been taken (e.g. Bonnell et al 1998, Myers 1998). The 
larger Av in regions forming massive stars imply that turbulence must be incorpo- 
rated into the models. Myers & Fuller (1992) suggested a "thermal-non-thermal" 
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(TNT) model, in which n(r) is represented by the sum of two power-laws and the 
term with p = 1 dominates outside the radius (ttnt) &t which turbulent motions 
dominate thermal motions. McLaughlin & Pudritz (1997) develop the theory of 
a logatropic sphere, which has a density distribution approximated by p = 1. 
Collapse in such a configuration leads to power laws in the collapsing region with 
similar form to those in the collapsing isothermal sphere, but with higher densi- 
ties and lower velocities. These ideas lead to accretion rates that increase with 
time and timescales for forming massive stars that are much weaker functions of 
the final stellar mass than is the case for the isothermal sphere (§ [4.5|) . Recent 
simulations of unmagnetized fragmentation that follow the interaction of clumps 
find that the mass spectrum of fragments steepens from a = 1.5 to a lognormal 
distribution of the objects likely to form stars (e.g. Klessen et al 1998). To avoid 
an excessive global star formation rate (§§) and distortion of the clump mass 
spectrum in the bulk of the cloud, this process must be confined to star-forming 
regions in clouds. 

5.4 Filaments, Clumps, Gradients, and Disks 

An important issue is whether the dense cores have overall density gradients or 
internal structures (clumps) that are likely to form individual stars and whether 
the mass distribution is like that of stars. Unfortunately, the terms "clumps" and 
"cores" have no standard usage; I will generally use "cores" to refer to regions 
that appear in maps of high-excitation lines and "clumps" to mean structures 
inside cores, except where a different usage is well established (e.g. "hot cores"). 
Myers (1998) has suggested the term "kernels" to describe clumps within cores. 
Cores themselves are usually embedded in structures traced by lines with lower 
critical density, and these are also called clumps by those who map in these 
lines. Many of these low-density structures are quite filamentary in appearance: 
examples include the 13 CO maps of L1641 (Orion B) of Bally et al (1987). In 
some cases, this filamentary structure is seen on smaller scales in tracers of high 
density or column density (e.g. Johnstone &; Bally 1999, Figure 2). 

The clumpy structure of molecular clouds measured in low-excitation lines 
suggests that dense cores will be clumpy as well. Suggestions of clumpy struc- 
ture came from early work comparing densities derived from excitation analysis 
in different tracers, but smooth density gradients provided an alternative (e.g. 
Evans 1980). Multitransition studies of three cores in CS (Snell et al 1984), C 34 S 
(Mundy et al 1986) and H2CO (Mundy et al 1987) found no evidence for overall 
density gradients; the same high densities were derived over the face of the core, 
while the strength of the emission varied substantially. This was explained in a 
clumpy model with clump filling factors of the dense gas /„ ~ 0.03 to 0.3, based 
on a comparison of My with M n (Snell et al 1984). This comparison forms the 
basis for most claims of unresolved clumps. 

Observations with higher resolution support the idea of clumps postulated by 
Snell et al (1984). For example, Stutzki and Glisten (1990) deconvolved 179 
clumps from a map of C 18 J = 2 — > 1 emission near M17. Because of overlap, 
far fewer clumps are apparent to the eye; assumptions about the clump shape 
and structure may affect the deconvolution. Maps of the same source in several 
CS and C 34 S lines (Wang et al 1993) could be reproduced with the clump catalog 
of Stutzki & Glisten, but only with densities about 5 times higher than they 
found. Thus the clumps themselves must have structure, either a continuation of 
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dumpiness or smooth gradients. Since the inferred clumps are now similar in size 
to the cores forming low mass stars, a natural question is whether massive cores 
are fragmented into many clumps which can be modeled as if they were isolated 
cores. In favor of this view, Stutzki and Giisten noted that the Jeans length was 
similar to the size of their clumps. 

A significant constraint on this picture is provided by the close confinement of 
the clumps; unlike the picture of isolated core formation, the sphere of influence 
of each clump will be limited by its neighbors. A striking example is provided 
by the dust continuum maps of the p Ophiuchi cloud (Figure 1), our nearest 
example of cluster formation, albeit with no very massive stars. Within about 
six cores of size 0.2 pc, Motte et al (1998) find about 100 structures with sizes of 
1000-4000 AU. They deduce a fragmentation scale of 6000 AU, five times smaller 
than isolated cores in Taurus. Thus the "feeding zone" of an individual clump is 
considerably less and the evolution must be more dynamic, with frequent clump- 
clump interactions, than is the case for isolated star formation. This picture 
probably applies even more strongly to the more massive cores. In p Ophiuchi, 
the clump mass spectrum above 0.5 M Q steepens to a = 2.5, close to the value for 
stars (Motte et al 1998), in agreement with predictions of Klessen et al (1998). A 
similar result (a = 2.1) is found in Serpens, using millimeter interferometry (Testi 
& Sargent 1998). While more such studies are needed, these results are suggesting 
that dust continuum maps do trace structures that are likely precursors of stars, 
opening the study of the origin of the IMF to direct observational study. 

Some of the less massive, relatively isolated cores, such as NGC2071, S140, 
and GL2591, have been modeled with smooth density and temperature gradients 
(Zhou et al 1991, 1994b, Carr et al 1995, van der Tak et al 1999). Models with 
gradients can match the relative strengths of a range of transitions with different 
excitation requirements, improving on homogeneous models. Zhou et al (1994) 
summarized attempts to deduce gradients and found preliminary evidence that, 
as core mass increases, the tendency is first toward smaller values of p. The most 
massive cores showed little evidence for any overall gradient and more tendency 
toward clumpy substructure. This trend needs further testing, but it is sensible 
if more massive cores form clusters. Lada et al (1997) found that the L1630 
cores forming rich embedded clusters with high efficiency tended to have larger 
masses of dense (n > 10 5 cm" 3 ) gas, but a lower volume filling factor of such gas, 
indicating more fragmentation. 

However, the line profiles of optically thick lines predicted by models with gra- 
dients are usually self-absorbed, while the observations rarely show this feature 
in massive cores. Clumps within the overall gradients are a likely solution. The 
current state of the art in modeling line profiles in massive cores is the work of 
Juvela (1997, 1998), who has constructed clumpy clouds from both structure tree 
and fractal models, performed 3-D radiative transport and excitation, and com- 
pared the model line profiles to observations of multiple CS and C 34 S transitions 
in massive cores. He finds that the clumpy models match the line profiles much 
better than non-clumpy models, especially if macroturbulence dominates micro- 
turbulence. Structure trees (Houlahan Sz Scalo 1992) match the data better than 
fractal models, but overall density and/or temperature gradients with p + q ~ 2 
are needed in addition to clumps. 

The study of gradients versus clumps in regions forming intermediate mass 
stars could help to determine whether conditions change qualitatively for star 
formation above some particular mass and how this change is related to the 
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outcome. Using near-infrared observations of regions with Herbig Ae/Be stars, 
Testi et al (1997) found that the cluster mode of star formation becomes dominant 
when the most massive star has a spectral type earlier than B7. Studies of the far- 
infrared emission from dust remaining in envelopes around Herbig Ae/Be stars 
found values of p ranging from 0.5 to 2 (Natta et al 1993). Maps of dust continuum 
emission illustrate the difficulties: the emission may not peak on the visible star, 
but on nearby, more embedded objects (Henning et al 1998, Di Francesco et al 
1998). Detailed models of several suitable sources yield p = 0.75 to 1.5 (Henning 
et al 1998, Colome et al 1996). Further work is needed to determine whether a 
change in physical conditions can be tied to the change to cluster mode. 

Many Herbig Ae stars have direct evidence for disks from interferometric stud- 
ies of dust emission (Mannings & Sargent 1997), though fewer Herbig Be stars 
have such direct evidence (Di Francesco et al 1997). For a review, see Natta et 
al (2000). Disks may be more common during more embedded phases of B star 
formation; Shepherd and Kurtz (1999) have found a large (1000 AU) disk around 
an embedded B2 star. The statistics of UC HII regions may provide indirect 
evidence for disks around more massive stars. Since such regions should expand 
rapidly unless confined, the large number of such regions posed a puzzle (Wood 
&z Churchwell 1989). Photoevaporating disks have been suggested as a solution 
(Hollenbach et al 1994). Such disks have also been used to explain very broad 
recombination lines (Jaffe & Martin-Pintado 1999). Kinematic evidence for disks 
will be discussed in § |5.5| . 

A particular group of clumps (or cores) deserve special mention: hot cores (e.g. 
Ohishi 1997). First identified in the Orion cloud (Genzel & Stutzki 1989), about 
20 are now known (see Kurtz et al 2000). They are small regions (I ~ 0.1 pc), 
characterized by Tk > 100 K, n > 10 7 cm -3 , and rich spectra, probably reflecting 
enhanced abundances, as well as high excitation (van Dishoeck & Blake 1998). 
Theoretical issues have been reviewed by Millar (1997) and Kaufman et al (1998), 
who argue that they are likely to be heated internally, but they often lack radio 
continuum emission. Since dynamical timescales for gas at such densities are 
short, these may plausibly be precursors to the UC HII regions. 

The evidence for flatter density distribution in regions of intermediate mass 
support the relevance of models like the TNT or logatropic sphere models in 
massive regions (§ |5.3|) , but it will be important to study this trend with the same 
methods now being applied to regions forming low mass stars, with due regard for 
the greater distance to most regions forming massive stars. The increasingly frag- 
mented structure in more massive cores and the increased frequency of clusters 
above a certain mass are consistent with a switch to a qualitatively different mode 
of star formation, for which different theories are needed. Finally, the common 
appearance of filaments may support a continuing role for turbulence in dense 
regions, since simulations of turbulence often produce filamentary structure (e.g. 
Scalo et al 1998). 

5. 5 Kinematics 

Lada et al (1991) found only a weak correlation between linewidth and size, dis- 
appearing entirely for a different clump definition, in the L1630 cores (Goodman 
Type 2 relation, see § |4.3| ) . Caselli & Myers (1995) also found that the Type 1 
linewidth-size relation (with only non-thermal motions included, Au/vt oc W) 
is flatter in massive cloud cores (7 = 0.21 ± 0.03) than in low mass cores (7 = 
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0.53 ± 0.07); in addition, the correlation is poor (correlation coefficient of 0.56) 
though the correlation is better for individual cores (Type 3 relations). They also 
noted that non-thermal (turbulent) motions are much more dominant in more 
massive cores and find good agreement with predictions of the TNT model. The 
"massive" cores in the Caselli & Myers study are mostly the cores in Orion with 
masses between 10 and 100 M . The much more massive ((My) = 3800 M ) 
cores studied by Plume et al (1997) exhibit no statistically significant linewidth- 
size relation at all (correlation coefficient is 0.26) and the linewidths are system- 
atically higher (by factors of 4-5) for a given size than would be predicted by the 
relationships derived for low and intermediate mass cores (Caselli & Myers 1995). 
In addition, the densities of these cores exceed by factors of 100 the predictions 
of density-size relations found for less massive cores (Myers 1985). The regions 
forming truly massive stars are much more dynamic, as well as much denser, 
than would be expected from scaling relations found in less massive regions. The 
typical linewidth in massive cores is 6-8 km s -1 , corresponding to a 1-D velocity 
dispersion of 2.5-3.4 km s _1 , similar to that of the stars in the Orion Nebula 
Cluster (Hillenbrand & Hartmann 1998). Larson (1981) noted that regions of 
massive star formation, like Orion and M17, did not follow his original linewidth- 
size relation, suggesting that gravitational contraction would decrease size while 
keeping Av roughly constant or increasing it. 

Searching for collapse in massive cores is complicated by the turbulent, clumpy 
structure, with many possible centers of collapse, outflow, etc. A collapse signa- 
ture may indicate an overall collapse of the core, with accompanying fragmenta- 
tion. In fact, self-absorbed line profiles from regions forming massive stars are 
rather rare (Plume et al 1997). A possible collapse signature has been seen in CS 
emission toward NGC 2264 IRS, a Class I source with L ~ 2000 L (Wolf-Chase 
k, Gregersen 1997). If an HII region lies at the center of a collapsing core, ab- 
sorption lines should trace only the gas in front and should be redshifted relative 
to the emission lines. Failure to see this effect, comparing H2CO absorption to 
CO emission, supported an early argument against the idea that all clouds were 
collapsing (Zuckerman & Evans 1974). A more recent application of this tech- 
nique to dense cores, using high-excitation lines of NH3, showed no preference for 
red-shifted absorption (Olmi et al 1993) overall, but a few dense cores do show 
this kind of effect. These sources include W49 (Welch et al 1988, Dickel & Auer 
1994), G10.6-0.4 (Keto et al 1988), and W51 (Zhang & Ho 1997, Zhang et al 
1998a). 

Dickel & Auer (1994) tested different collapse scenarios against observations 
of HCO + and favored free-fall collapse, with n cx r~ 15 and v oc r~°- 5 throughout 
W49A North; they noted that more complex motions are present on small scales. 
Keto et al (1988) used NH3 observations with 0.3" resolution to separate infall 
from rotational motions toward the UC HII region, G10.6-0.4. Zhang & Ho 
(1997) used NH3 absorption and Zhang et al (1998a) added CS J = 3 — ► 2 and 
CH3CN observations to identify collapse onto two UC HII regions, W51e2 and 
W51e8, inferring infall velocities of about 3.5 km s _1 on scales of 0.06 pc. Young 
et al (1998) have tested various collapse models against the data on W51e2 and 
favor a nearly constant collapse velocity (y ~ 5 km s _1 ) and n(r) oc r~ 2 . Mass 
infall rates of about 6 x 10" 3 (Zhang et al 1998a) to 5 x 10" 2 (Young et al 1998) 
M yr _1 were inferred for W51e2. Similar results were found in G10.6-0.4 (Keto 
et al 1988), and even more extreme mass infall rates (10~ 2 to 1 M yr _1 ) have 
been suggested for W49A, a distant source with enormous mass (~ 1O 6 M ) and 
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luminosity {L ~ 10 7 Lq) (Welch et al 1988). These high infall rates may facilitate 
formation of very massive stars (§ |5.1[ ) and help confine UC HII regions (Walmsley 
1995). 

Large transverse velocity gradients have been seen in some hot cores, including 
G10.6-0.4 (Keto et al 1988), W51 (Zhang k Ho 1997), and 77L4S20126 + 4104 
(Cesaroni et al 1997, Zhang et al 1998b). For example, gradients reach 80 
km s" 1 pc" 1 in the NH 3 (J, K) = (4, 4) line and 400 km s _1 pc _1 in the CH 3 CN 
J = 6 -> 5 line toward G29.96-0.02 and G31.41+0.31 (Cesaroni et al 1994b, 
1998). Cesaroni et al (1998) interpret the gradients in terms of rotating disks 
extending to about 10 4 AU. Using the dust continuum emission at 3 mm, they 
deduce disk masses up to 4200 M Q , in the case of G31.41+0.31. 

5. 6 Implications for Larger Scales 

Since most stars form in the massive cores discussed in this section (Elmegreen 
1985), they are most relevant to issues of star formation on a galactic scale. If 
the luminosity is used to trace star formation rate (e.g. Rowan-Robinson et al 
1997, Kennicutt 1998), the star formation rate per unit mass is proportional to 
L/M. Considering clouds as a whole, {L/M) ~ 4 in solar units (e.g. Mooney & 
Solomon 1988), with a spread exceeding a factor of 10 2 (Evans 1991). Using CS 
J = 5 — ► 4 emission to measure M in the dense cores, Plume et al (1997) found 
(L/M) = 190 with a spread of a factor of 15. The star formation rate per unit 
mass is much higher and less variable if one avoids confusion by the sterile gas. 
The average L/M seen in dense cores in our Galaxy is similar to the highest value 
seen in ultra-luminous infrared galaxies, where M is measured by CO (Sanders et 
al 1991, Sanders & Mirabel 1998). The most dramatic starburst galaxies behave 
as if their entire interstellar medium has conditions like those in the most active, 
massive dense cores in our Galaxy. Some studies (e.g. Mauersberger & Henkel 
1989) indeed found strongly enhanced CS J = 5 — > 4 emission from starburst 
galaxies. It will be interesting to observe high-J CS lines in the most luminous 
galaxies to compare to the conditions in massive cores in our Galaxy. Perhaps 
the large scatter in L/M seen in galaxies will be reduced if CS, rather than CO, 
is used as a measure of the gas, ultimately leading to a better understanding of 
what controls galactic star formation rates (see Kennicutt 1998). 

5. 7 Summary of Clustered Star Formation 

The cloud mass distribution found for lower mass objects continues to massive 
clouds, but less is known about the distribution for dense cores. Very massive 
cores clearly exist, with sufficient mass (M > 10 3 M Q ) to make the most mas- 
sive clusters and associations. These cores are denser and much more dynamic 
than cores involved in isolated star formation, with typical n ~ 10 6 cm -3 and 
linewidths about 4-5 times larger than predicted from the linewidth-size rela- 
tion. Pressures in massive cores (both thermal and turbulent) are substantially 
higher than in lower mass cores. The densities match those needed to form the 
densest clusters and the most massive stars by coalescence. There are some cores 
with evidence of overall collapse, but most do not show a clear pattern. There is 
some evidence that more massive regions have flatter density profiles, and that 
fragmentation increases with mass, but more studies are needed. High resolution 
studies of nearby regions of cluster formation are finding many clumps, limiting 
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the feeding zone of a particular star-forming event to I ~ 6000 AU, much smaller 
than the reservoirs available in the isolated mode. In some cases, the clump mass 
distribution approaches the slope of the IMF, suggesting that the units of star 
formation have been identified. Studies of intermediate mass stars indicate that 
a transition to clustered mode occurs at least by a spectral type of B7. 

6 CONCLUSIONS AND FUTURE PROSPECTS 

The probes of physical conditions have been developed and are now fairly well 
understood. Some physical conditions have been hard to measure, such as the 
magnetic field strength, or hard to understand, notably the kinematics. Star- 
forming structures, or cores, within primarily sterile molecular clouds can be 
identified by thresholds in column density or density. Stars form in distinct 
modes, isolated and clustered, with massive stars forming almost exclusively in 
a clustered mode. The limited number of measurements of magnetic field leave 
open the question of whether cores are subcritical or supercritical, and whether 
this differs between the isolated and clustered mode. Cores involved in isolated 
star formation may be distinguished from their surroundings by a decrease in 
turbulence to subsonic levels, but clustered star formation occurs in regions of 
enhanced turbulence and higher density, compared to isolated star formation. 

Evolutionary scenarios and detailed theories exist for the isolated mode. The 
theories assume cores with extended, power-law density gradients, leading to the 
mass accretion rate as the fundamental parameter. Detailed tests of these ideas 
are providing overall support for the picture, but also raising questions about 
the detailed models. Notably, kinematic evidence of gravitational collapse has 
finally been identified in a few cases. The roles of turbulence, the magnetic field, 
and rotation must be understood, and the factors that bifurcate the process into 
single or multiple star formation must be identified. Prospects for the future 
include a less biased census for cores in early stages and improved information 
on density gradients, both facilitated by the appearance of cameras at millimeter 
and submillimeter wavelengths on large telescopes. Antenna arrays operating at 
these wavelengths will provide more detailed information on the transition region 
between envelope and disk and study early disk evolution and binary fraction. 
Future, larger, antenna arrays will probe disk structure to scales of a few AU. 
Finally, a closer coupling of physical and chemical studies with theoretical models 
will provide more pointed tests of theory. 

Theories and evolutionary scenarios are less developed for the clustered mode, 
and our understanding of the transition between isolated and clustered mode is 
still primitive. Current knowledge suggests that more massive cores have flatter 
density distributions and greater tendency to show substructure. At some point, 
the substructure dominates and multiple centers of collapse develop. With re- 
stricted feeding zones, the mass accretion rate gives way to the mass of clump 
as the controlling parameter, and some studies of clump mass spectra suggest 
that the stellar IMF is emerging in the clump mass distribution. The most mas- 
sive stars form in very turbulent regions of very high density. The masses and 
densities are sufficient to form the most massive clusters and to explain the high 
stellar density at the centers of young clusters. They are also high enough to 
match the needs of coalescence theories for the formation of the most massive 
stars. As with isolated regions, more and better measurements of magnetic fields 
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are needed, along with a less biased census, particularly for cool cores that might 
represent earlier stages. Larger antenna arrays will be able to separate clumps 
in distant cores and determine mass distributions for comparison to the IMF. 
Larger airborne telescopes will provide complementary information on luminos- 
ity sources in crowded regions. Observations with high spatial resolution and 
sensitivity in the mid-infrared will provide clearer pictures of the deeply em- 
bedded populations, and mid-infrared spectroscopy with high spectral resolution 
could trace kinematics close to the forming star. Deeper understanding of clus- 
tered star formation in our Galaxy will provide a foundation for understanding 
the origin and evolution of galaxies. 
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Figure 1: (NOTE: this figure is not supplied here. It can be seen in the original 
paper (ref. below) and it exists on astro-ph (9903284) as a gzipped postscript file.) 
Map of 1.3 mm dust continuum emission covering about 2 pc in the p Ophiuchi 
cloud, with a resolution of 15" (Motte et al 1998). The emission is proportional 
to column density, and the lowest contour corresponds to N ~ 1.4 x 10 22 cm" 2 
or Ay ~ 14 mag. Contour levels increase in spacing at higher levels (see Motte 
et al 1998 for details). The highest contour implies Ay ~ 780 mag. 



Figure 2: (NOTE: this figure is not supplied here. It can be seen in the original 
paper (ref. below) in grayscale.) The 850 fim emission and 450 fj,m to 850 /im 
spectral index distribution, covering about 7 pc with a resolution of 14" at the 
northern end of the Orion A (L1641) molecular cloud (Johnstone & Bally 1999). 
The Orion Nebula is located directly in front of the strong emission at the center 
of the image. (Left) The 850 /im image showing the observed flux from —0.1 to 
2 Jy/beam with a linear transfer function. The highest flux level corresponds 
roughly to Ay ~ 320 mag if Tk = 20 K. (Center) The 850 /im image showing 
the observed flux from 100 mJy/beam to 20 Jy/beam with a logarithmic transfer 
function. (Right) The 450 /jm to 850 /im spectral index in the range 2 to 6 with 
a linear transfer function. 



